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ABSTRACT 


A large-signal theory of oscillator noise is presented, which is 
applicable to oscillators employing active elements that can be described 
by a time-dependent differential equation with a Langevin noise source. 

The method uses a modified quasistationary perturbation theory with an RF 
amplitude-dependent perturbing noise voltage (current) source, which has 
been mathematically expressed in terms of the "initial" (pre-oscillation) 
amplitude- and phase noise modulation rates. The usual requirement that 
this noise voltage (current) be a member of a normal random process, which 
results in equality of the two initial power spectra, is not invoked here. 
The RF amplitude dependence of the open-circuit perturbing noise voltage 
(short-circuit noise current) has been obtained from the differential 
equation of the active element, under the condition of sinusoidal current 
(voltage) excitation at the frequency at which it would normally oscillate. 
Upon linearizing in stochastic quantities, expanding the deterministic terms 
by means of Fourier series and Fourier transforming the result, a system 

of linear equations in the frequency domain is obtained for the noise 
voltages (currents) at frequencies about multiples of the excitation 
frequency and at the baseband frequency. In this way, most of the conversion 
processes in the active element can be taken into account. The Fourier 
transforms and power spectra of the initial modulation rates have been 
determined in terms of the upper- and lower sideband noise voltages 
feacrenteh: they are used in the modified formulas of a quasistationary 
perturbation theory to obtain the AM-and FM noise spectra of the oscillator 


output. 
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The theory has been applied to IMPATT diode oscillators and the 
computed spectra have been compared with various experimental data. To 
supplement data already available in the literature, a sensitive noise- 
measuring system has been assembled, carefully tested and used for 
amplitude- and frequency noise measurements on IMPATT diodes. The 
experimental data have been found to be in good agreement with the theore- 
tical predictions. | 

Additional insight into IMPATT diode noise has been provided, since 
the above method made it possible to evaluate the individual (and generally 
different) dependence of the AM-and FM noise spectra on the RF output 
power of the oscillator as well as on other diode- and circuit parameters. 
As a result of the theoretical and experimental studies, guidelines have 
been obtained for determining the optimum parameters of an IMPATT diode 


oscillator, for minimum noise. 
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CHAPTER I 
INTRODUCTION 


The subject of oscillator noise, its theoretical description and 
experimental assessment is not only an interesting problem per se, but 
is also one of direct practical importance. The amount of oscillator- 
generated noise determines to an important extent the quality of a 
telecommunications system. The invention of a new electron device 
that is capable of operating as an oscillator has therefore generally 
been followed by noise studies of the active element itself and of the 
oscillator unit in which the latter is used. This has been especially 
true in the case of recent semiconductor devices, such as Gunn diodes, 
IMPATT diodes and BARITT diodes, all of which are capable of oscillating 
in the microwave frequency range. Specifically, IMPATT diodes with their 
strong nonlinearities spurred the interest in large-signal studies of 
oscillator noise. In this thesis, it has been attempted to develop a 
large-signal theory of oscillator noise that would take into account most 
of the conversion processes existing in high-level oscillators, while at 
the same time yielding a clearly defined procedure for actual. computation 
of the oscillator noise spectra, the results of which could be compared 


with measured data. 
1-1 Brief Review of Pertinent Oscillator Noise Theories 


The theory of noise in free-running oscillators was first established 
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by Bernstein’. From the studies that followed, the ones most often cited 
and actually used for oscillator noise assessment are those of Edson“, 


757) These theories provide relatively simple 


Mullen® and Kurokawa 
analytic expressions for the amplitude (AM)-and frequency (FM) noise 
spectra. 

In order to obtain these simple analytic expressions, the following 
basic assumptions needed to be made in the above studies: 1) the amplitude- 
and phase fluctuations are small, so that a quasilinear approximation is 
possible about the operating point of the active element, 2) the random 
amplitude and phase vary slowly with time (quasistationarity), 3) the 
perturbing noise voltage source used in the oscillator equivalent circuit 
is independent of the amplitude of oscillation in a small region near the 
operating point, and 4) the perturbing noise voltage is decomposed into 
two components at the frequency of oscillation, one in-phase and the other 
in quadrature with the oscillator periodic signal; the slowly varying 
random amplitudes of these components are assumed to be uncorrelated and 
to have equal power spectra. In addition, the influence of the oscillation 
frequency harmonics on the noise is neglected. 

A more rigorous and systematic study by Lax® has shown that, while 
the quasilinear approximation is valid for amplitude fluctuations (if this 
approximation is made in the correct variable), phase fluctuations should 
not be treated in this manner, since they may be quite large in self- 
sustained oscillators. Beside other results, Lax's paper contains 


comparisons of spectral linewidths obtained by different procedures. 


Hempstead and Lax! go further in this by exactly solving the appropriate 
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Fokker-Planck equation, which is probably the most rigorous approach to 
take when the fluctuations are so large that a quasilinear approximation 
is not warranted. 

Although the above considerations are of fundamental importance, 
they do not have a crucial effect in the practical application of a noise 
theory to typical oscillators in the radio-frequency and microwave range. 
This is because, even under large-signal conditions, the fluctuations in 
phase and in amplitude are rarely so large at these frequencies as to 
seriously invalidate the quasilinear approximation. While this approximation 
is incorrect in principle for phase fluctuations, its quantitative effect 
is manifested strongly mainly in the spectral linewidth, and much less so 
in the spectral density of frequency fluctuations; this quantity is pre- 
dominantly used for frequency noise characterization of the above class of 
oscillators. It has been estimated by this author, using the appropriate 
_ formulas of Lax®, that in a microwave oscillator at 10 GHz, a 3 dB 
difference in the spectral density of frequency noise caused by the quasi- 
linear approximation will typically occur at very low baseband frequencies, 
100 Hz - 1 kHz. This is a negligible portion of the baseband frequency 
region of interest. In other respects, Lax's theory® is not fundamentally 
different from the previously mentioned methods; it uses the assumption 
of quasistationarity and that of a spectrally pure oscillator. The 
perturbing noise voltage is formally dependent on a controlling parameter, 
associated with the amplitude of oscillation; understandably, since the 
theory is a general one, no details are given as to how this dependence 


arises in a specific case. 
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As a consequence of the lack of a complete and computationally 
applicable large-signal theory, the assessment of noise in typical 
oscillators used for telecommunications purposes has so far been based on 
Edson's and Kurokawa's results, in which the perturbing noise voltage is 
normally replaced by the noise voltage present across the active element 
open-circuited at the frequency of oscillation and dc-excited only. For 
devices having a weak dependence of the open-circuit noise voltage on the 
ac voltage and current that are present in the circuit under oscillatory 
conditions, this approach is satisfactory. 

With the commencement of noise studies on IMPATT diode oscillators, 
jt soon became apparent that some of the noise properties of oscillators 
employing such a highly nonlinear element will defy description by the 
existing theories. A rather extensive literature has developed, dealing 
with both small-signal and large-signal noise properties of IMPATT diodes. 
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After the fundamental studies of Tager’, Hines~, Haitz 


1] 


and Gummel and 


Blue on dc-excited IMPATT diodes, attention shifted to the phenomenon 


of up-conversion and its different influence on the amplitude-and frequency 


lz, and to the means of suppressing the 


13) 


spectra, first treated by Scherer 


frequency noise by injection locking (Kurokawa and Hines et al. or 


V4) Many other authors 


by high-Q cavity stabilization (Ashley and Searles 
studied various aspects of the IMPATT noise problem. Rather than present 
the extensive list of references here, we refer the reader to the review 
paper by Gupta!>, where most of the contributions in this field are 
briefly discussed and a list of references is given. However, some of 


these studies, considered by the present author to be of direct importance 


to the subject of this thesis, will now be briefly reviewed. 
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Thaler et cata? presented a quasistationary perturbation theory, 
based on Kurokawa's work, in which the influence of bias current © 
fluctuations on the diode impedance at the oscillation frequency is 
introduced. This permitted the expression of amplitude-and frequency 
noise in terms of primary- and modulation noise. Although this theory 
is a small-signal one, since the perturbing noise voltage is assumed to 
be independent of the amplitude of oscillation, it yields some interesting 
suggestions on minimizing the internal phase-to-amplitude conversion in 
the oscillator (or phase-to-amplitude coupling, as mentioned by Lax®). 

convert! 7 is probably the first to have pointed out that the 
perturbing noise voltage in IMPATT diode oscillators may be strongly 
dependent on the amplitude of oscillation and that its two components, 
in- and out-of-phase with the oscillator signal, may be cross-correlated. 
His derivations were carried out specifically for the Read diode. The 
disadvantage of this rather involved treatment is that it puts aside some 
of the circuit aspects; also, no computational results are available that 
would permit comparisons with measured data. The theory uses a previously 
published large-signal study of the Read diode by Delagebeaudeuf |®, 

Vlaardingerbroek !? modified the quasistationary perturbation theory 


of Kurokawa’ ?? 


by including the effect of down-conversion by demodulation. 
The resulting system of equations is the most complete one in the lit- 

erature as of this writing, and provides a sound basis for a large-signal 
generalization if the noise source voltages can be computed in dependence 
on the amplitude of oscillation. In such a generalization, the amplitude 


dependence of the noise source voltages and their cross spectra, normally 


neglected in the small-signal approach, would have to be considered. 
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A very important step in the analysis of the amplitude dependence of 


noise in IMPATTs is the study of Kuvas??. 


First, the statistics of the 
basic noise-generating process, namely that of random ionization in the 
avalanche region, are derived in dependence on the amplitude of the driving 
current at the frequency at which the device would normally oscillate. 
He concluded that, for the usual range of excitation, this influence is 
negligible and the spectrum density of this basic noise process depends 
only on the magnitude of the bias (dc) current, as first derived in the 
small-signal treatment of Tager®. Next, the transformation of this random 
process into its external form, the open-circuit noise voltage, is described 
mathematically under conditions of large-signal excitation. It is shown that 
the open-circuit noise voltage is strongly dependent on the magnitude of 
the periodic driving force, although the basic noise process in the avalanche 
region is little influenced by it. Alternatively, the noise measure and 
power-to-noise ratio are calculated in dependence on the oscillator RF power. 
All calculations are carried out for an infinite impedance in the bias circuit 
of the IMPATT diode, so that modulation noise caused by random noise current 
at low frequencies is excluded. It is suggested by the author that some 
information on the AM-and FM noise spectra can be obtained from his theory 
in conjunction with perturbation theories of oscillator noise, such as those 
in Refs.2, 4, and 5. 

This author wishes to emphasize that the use of noise measure gives 
correct results for AM-and FM noise only in the case where the slowly- 
varying amplitudes of the in- and out-of-phase components of the perturbing 


noise voltage have equal spectra. This is typical for low-level oscillators. 
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If this concept is used as suggested in Ref. 20, the increase in noise 
measure with the amplitude of oscillation will be reflected in AM-and 

FM noise spectra in a precisely identical manner. This would be in 
complete disagreement with experimental data, e.g., in Ref. 21, where 

the AM noise of IMPATT diodes has been found to change much less with 

the RF voltage across the diode than FM noise, in the region of high RF 
voltage. The explanation fede on the fact that AM-and FM noise spectra 
are uniquely determined by both the random amplitude and phase of the 
perturbing noise voltage sidebands above and below the frequency of 
oscillation. However, Kuvas' procedure computes only the spectrum density 
of the noise voltage at the oscillation frequency, i.e., at zero 
modulation frequency. With this approach, it is impossible to distinguish 
between the generally different RF amplitude dependence of the in- and 
out-of-phase components of the perturbing noise voltage, and that of their 
_ auto- and cross spectra. As a consequence, the AM-and FM noise spectra of 
the oscillator output cannot be correctly computed. To obtain separate 
expressions for AM-and FM noise spectra, a different approach must be 
employed right at the beginning of large-signal derivations, although some 
of Kuvas' ideas, such as application of Fourier transform and linearization 


in noise amplitudes can be successfully utilized. 


1-2 Purpose of this Study 


A critical assessment of the available oscillator noise theories 


Shows the need for a complete large-signal theory that would bridge the 
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existing gaps and provide a tool for a sufficiently accurate calculation 
of noise in typical radio-frequency oscillators. It is the purpose of 
this study to develop a theory of oscillator noise that will consistently 
take into account the joint influence of the amplitude of oscillation and 
that of the higher harmonics, of the active element parameters and that 
of external circuitry on the output amplitude-and frequency noise. 

The theory is intendedmne be applied to IMPATT diode oscillators 
and the computed spectra will be compared with experimental data. Results 
of the measurements of other authors will be supplemented by the 
experimental data of the present author which were obtained on IMPATT 
diode oscillators by means of a sensitive noise-measuring system, which 
was assembled specifically for this project. It is hoped that this work 
will help elucidate some of the noise phenomena in IMPATT diodes which are 
already known, and also provide a more general insight into IMPATT 

oscillator noise. Finally, the quantitative comparisons of the computed 
| and measured spectra should yield answers as to the viability of the 


proposed theory. 


1-3 Organization of this Thesis 


The main part of the original material is contained in Chapters 
II to VI and some auxiliary work is presented in Appendices A to C. 
Chapter 1] is devoted to the discussion- and large-signal modification of 
the existing quasistationary perturbation theories of oscillator noise. 
The concept of pre-oscillation modulation rates is introduced in this 


chapter. In Chapter III, a general large-signal procedure is developed 
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for the computation of the open-circuit noise voltage and of the pre- 
oscillation modulation rates for the active element. Chapter IV contains 
the explicit application of the above general procedures to IMPATT diodes, 
complete with the plots and discussion of the computed amplitude-and phase 
(frequency) noise spectra. Some auxiliary derivations and formulas 
supporting the developments of Chapter IV are given in Appendices A to C. 

A comparison of the Mente spectra with experimental data is 
carried out in Chapter V; the noise-measuring system used is also 
described here. 


A summary and conclusions are presented in Chapter VI. 
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CHAPTER II 
QUASISTATIONARY PERTURBATION THEORY OF NOISE IN FREE-RUNNING 
OSCILLATORS: MODIFICATION FOR LARGE-SIGNAL CONDITIONS 


One of the aims of the large-signal treatment of oscillator noise 
proposed in this thesis is the removal of certain restrictions imposed 
on the perturbing noise source in the theories cited in Chapter I. A 
Suitable "small-signal" perturbation theory can then be modified 
accordingly so that the resulting expressions for oscillator output 
noise can accept, as an input, quantities characterizing the amplitude 
dependence of the perturbing noise source. The above concept is 
elaborated on in this Chapter, utilizing the "state of the art" in 
quasistationary perturbation theory as represented by Thaler et ae 


and Viaardingerbroek!9, whose contributions are based on the fundamental 


paper of Kurokawa’. 


2-1 Principles of the Perturbation Method in Oscillators 


2-1-1 Oscillator Model with Perturbing Noise Source 


For purposes of noise analysis in general, the equivalent circuit 
of an BeciLlaroe should include both the high-frequency (RF) part as 
well as the low-frequency (baseband) part, through which the bias voltage 
and current are also supplied to the active element. Such a combined 


circuit for oscillator noise analysis has been introduced only recently 
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19° This circuit is shown, in a somewhat modified 


by Vlaardingerbroek 
form, in Fig. 2-la. The impedance of the active element in the RF 
circuit is,in general ,dependent on the frequency w(t) and the amplitude 
A‘ (t) of high-frequency oscillations as well as on the bias current 
I(t). In the low-frequency circuit, the impedance of the active element 
depends also on the frequency 2. The perturbing voltage source v(t) is 
taken to represent the noise process in the active element. For the 
present purpose this source may be separated into two sources, one for 
the baseband frequency range and, the other for the oscillation frequency 
range. These two voltages may, in general ,be correlated (see Sec. 2-2-2). 
For completeness, we note that the active element can also be 
represented by the dual of Fig. 2-la, an admittance in parallel with a 
current noise source; the corresponding equivalent circuit is in Fig. 2-1b. 
The behaviour of these two circuits is Similar as far as the frequency 
_ region around the fundamental frequency of oscillation is concerned; the 
main difference is that, for oscillators with a sufficient quality factor 
Q, the RF circuit in Fig. 2-la contains voltage components at higher 
harmonics but only one current component (at the fundamental frequency) 
whereas, in Fig. 2-1b, the situation is reversed. This difference is of 
no importance in a small-signal noise analysis, where the perturbing 
noise source is independent of the magnitude of the deterministic current 
or voltage at the active element. However, as will become apparent in 
Chapter III, in a large-signal calculation of the perturbing noise voltage 
or current the differences in harmonic content may play a significant 
role. It is, therefore, important to choose the correct equivalent 


circuit for the oscillator to be investigated. In this thesis we shall 
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Fig. 2-la OSCILLATOR EQUIVALENT CIRCUIT WITH NOISE VOLTAGE SOURCE 
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Fig. 2-lb OSCILLATOR EQUIVALENT CIRCUIT WITH NOISE CURRENT SOURCE 
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give specific relations only for the circuit in Fig. 2-la, as this is 


the one of importance for our intended application. 


2-1-2 Representation of the Oscillator Output Signal and of the 
Perturbing Noise Source; Output Noise Modulation Rates 


The current in the RF loop of the oscillator circuit in Fig. 2-la 
is assumed to be concentrated in frequency around a fixed frequency of 


oscillation Wo and it can therefore be expressed as follows: 


I(t) = A'(t)cosa(t) = [At + aA'(t)]cos[w,t + na(t)] 


i] 


atl + M(t)]eos[w t + as(t)] (2.1) 


=< 

— 
c+ 

sd 
| 


: aa'(t)/Al (2.2) 


_ The output voltage, i.e., the voltage across the real part of the load 
impedance, is then given by 

uae 

= AR, (w i (223) 


| ae Vf + M(t) ]cos[w t + ao(t)] 3; V 


Re Cc Cc 
The frequency We is commonly called "carrier frequency". For the 
purpose of a noise analysis, no intentional deterministic modulation is 
assumed; both M(t) and Ae(t) are stochastic quantities representing the 
random fluctuations of the amplitude and phase of the oscillator signal. 
These fluctuations are caused by the presence of the noise sources in 


the oscillator circuit. We shall call M(t) and Ae(t) the output amplitude- 


and phase noise modulation rates, respectively. The power spectra of 


“| ef , yi 
fi ‘i oe 
ads Posbae ian? suo, ram A 2.9 ib B13 
ogden » po ith ube) a. 92" oie buh rc 0. PED aod gate : 
at O14 nt dias fo radi Eaten, to doo oa mo pm 
Fo yYonsupay? bext? a! paver arieupey oe pstgreanion ado 
:2wolto® 26 benaaairip od srohenstly 09 
[(3 env a, asia ha" ws ath = on 
(ay © 4 A = Wea 7 » Eeehoe a 
2 a a as 
alt 103 .*yonaupen?) ‘bisa ‘bel igo walle 3! | alleged it 
‘ey 
ef cohtet soon 914 ethane set “Uaholtnaba “on 2a lene sefon's YO "S20m 
anit: patinsesrqat’ eatt iansup atgeplionte ‘iB Cayee bine (sm tied coamvet 
sTaiipke, rodeb M20, ado oH Hs at ‘ott to eaten asl a 3 


~f ‘ fit zoriuoz 2rGn ‘ont Tol sonszbing att ad bpetisc 278 
(eee ara ait nee Gin fs (torte sie. .thadto: 


, ay ar ; 


14 


these quantities are extensively used for the noise characterization 

of oscillators and are referred to as "double-sideband amplitude (AM) 
noise-to-carrier ratio" and "double-sideband phase (PM) noise-to-carrier 
ratio’, both evaluated in a certain specified frequency band, e.g., 

100 Hz (located at a specified "distance" from weds 


In the important quasistationary limit, i.e., for M(t) and Ao(t) 


varying slowly with time, the output signal can be expressed in the form 


YR, vans M(t) Jcos{[w, Sehaey ley, ts (2.4) 


where Aw(t) is the random frequency deviation due to noise. The square 
root of its power spectrum is also extensively used and is referred to 
as the rms frequency deviation in a certain specified frequency band. 


Note that, in general, 


aw(t) = aLde(t)] ees) 


or in the frequency domain, 


Aw(2) = jao(n) , (2.6) 
so that the power spectra S yy 6) and S4q(%) of Aw(t) and Ao(t), 
respectively, are related by 

SP) Cee 0 S59 (2) ‘ (2.7) 


For completeness, we note that relations similar to Eqs. (2.1) - (2.4) 


can be written for the circuit in Fig. 2-1b. 
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In correspondence with the mathematical form used for the output 
signal, the perturbing noise voltage (current) in the RF region is 
conveniently expressed by means of the well-known narrowband represent- 


ation !® 


WwW = cy * 
v(t) v.(t)coswt v.(t)sinw t 


iH] 


R,{[v.(t) + jv. (t)Jexpjw t} ; (2.8) 


where Mek ty) and v,(t) are small stochastic quantities slowly varying in 
time. A detailed examination of the properties of these quantities will 


be carried out in Sec. 2-2. 


2-1-3 Quasistationary Perturbation Equations 
16 


The conditions of quasistationarity can be stated as follows -: 


clan (e)1. Mw ; & du(t) << w? : ee.) 
This means that All (t) and Aw(t) are assumed to vary only slowly with 
time. Such a restriction can be imposed if the region of interest is 
confined to low modulation frequencies in the power spectra of AA‘ (t) 
and Aw(t). 

This quasistationary approximation makes it possible in the RF 
region to express the voltage across (current through) the active element 
as the product of current (voltage) and active element impedance 


(admittance) at the instantaneous amplitude and frequency and also bias 


current. Similarly, the voltage across (current through) the load 
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impedance (admittance) is given as the product of current (voltage) 
and load impedance (admittance) at the instantaneous frequency |°. As 
a consequence, the voltage balance equation in the RF loop of Fig. 2-la 


reduces from a differential equation to the simple relation 


(A (t)expjlu, + au(t)]t)-(ZeCA" (t) ,1,(t) so(t)] + Z,[o(t) I} 
+ [v (t) + jv, (t) ]expjut S0e 0 8 (2.10) 


Now the deviations aA'(t), Al, (t) and Aw(t) from the stable 
operating point given by Al ln and w, are assumed small and the impedances 
in Eq. (2.10) are expanded into Taylor series, stopping with the linear 


terms. The following general perturbation equation results: 


i ; yes . a7” 
(TA. + AA (t) JexpjLu, + dw (t)]t}-{—# -AA'(t) a wg. % AI, (t) 
Cc bure 
a i : rat 
+ TPL, to dn, w(t)? + [velt) + Ivg(t)lexpiu,t = 0 


(211) 


Using the condition for steady-state oscillations 

Wypl a 

Za(Acs Ips wo) + Z, (w.) SO (22) 
neglecting the higher-order terms (products of small stochastic quanti- 
ties) in Eq. (2.11) and separating Eq. (2.11) into real and imaginary 


parts yields two perturbation equations 


middie aids worl (3a bas (19,14 sighed eee ‘a 
, @sansbaqmt’ ord brs Thome bemuzes” ox5 5 DAB gl Sow eave ott 


“eon! edd AdtWw pabagote «29482! xol ys} aga ‘tne 9 cons 
:2aTu2ey nottsups. notiedtussa revere 


4p) ghh 7s + (de. 
Ya ; . fc ha j 
i ) -* aa 


- = Hiasup si sesdoode: ere 12 bo 
‘ i  Nveatoent bne ts on (11 


17, 


oR ary ara dR v(t 
Ai AA (t) + i Al, (t) * Le + Gp JAw(t) fae (2.18) 
fe) res A 
Cc C 
So ax oye, med X v_(t) 
—2 on'(t) + Sal, (t) + [52+ Jaw(t) -- . (2.14) 
oA 0 aCe A. 


In these equations, the following notation has been introduced for 


convenience: 
ORG 9R® aR aR” 
aks) Beige!) 0h 2 ed: (2.15) 
All 4 tartan = | Shee 
A =A re I,=I 
Cc DieO 
etc. 


Eqs. (2.13) and (2.14), transformed into the frequency domain, have the 


form 
aR® aR” aR” dR, v_(2) 
ek Ani add -nlein ie—a ott (Q)aret ee (2.16) 
al b ow dw 1 
oA ) C Cc A. 
ce a a Bead i v_ (2) 
gE Hid (ded, aoe ott, shear Hee Dn “ae (eds ayia ee. ae ee) 
ol b ow dw 1 
3A 0) C ¢ A. 
where Vo(2)s Vo (2)5 wheel ek are the Fourier transforms of the correspond- 


ing quantities in the time domain, with the implicit assumption that the 
latter are truncated on a finite time interval (-T,T) and that they 


satisfy all the other requirements for the existence of the Fourier 


transform<. 


To complete the description of the oscillator circuit in Fig. 2-la, 


a bias circuit equation is introduced! ?: 
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y(ay(zg(a’, Ips, 2) +Z(0)} + v(a)+Vi = 0 . (2.18) 


d 
The meaning of the quantities contained in Eq. (2.18) is apparent from 
Fig. 2-la. Again, the impedance of the active element is expanded into 
a Taylor series around the operating point and, after removing the dc 


part, the following perturbation equation is obtained: 


az" a7" 
1g \al Oiea ae 
I, 3A AA (2) “uy [ZA Ths soe Q) + aT, aa Z(2) Jal, (2) 
c 
+ i am Aw(2) = - v(2) (2.19) 


The first and third terms in this equation represent the effect of 
demodulation!?. 

It may be useful at this point, and for future reference also, to 
depict the interaction between the individual quantities in Eqs. (2.16), 
(2.17) and (2.19) graphically. This is illustrated in the interaction 
diagram", Fig. 2-2, which shows the links between physical causes and 
effects. Note that, if modulation is nonexistent (224/21, = 0), AA! (a) 


and Aw(@) are independent of v(2). Similarly, v.(2) or V.(@) ean 


S 
influence AT, (2) only via demodulation. 

If the equivalent circuit in Fig. 2-lb were used, the system 
corresponding to Eqs. (2.16), (2.17) and (2.19) would be completely 
analogous. The noise voltage sources would be replaced by noise current 


sources, the bias current fluctuations AI, (t) by voltage fluctuations 


AV, (t) and all impedances by admittances. 
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Mathematically, Eqs. (2.16), (2.17) and (2.19) form a linear system 
that can be solved for AA' (a) ,Aw(2) and AT, (2) in terms of the perturbing 
voltages V (2). v. (2) and v(2). This has been carried out by Thaler et 


al gto 


and Viaardingerbroek!9, using different simplifying assumptions; a 
brief review of this is given in Sec. 2-3. Prior to this, the spectral 


properties of the perturbing noise source must be investigated. 


2-2 Properties of the Perturbing Noise Source 


2-2-1 Introduction of the Initial (pre-oscillation) Noise 
Modulation Rates 

The form of Eq. (2.8) indicates that the perturbing noise voltage 
could be expressed by means of some equivalent modulation rates or 
indices characterizing the initial (pre-oscillation) amplitude-and phase 

_ fluctuations. This would also formally complement the introduction of 

noise modulation rates in the output noise characterization, Eqs. (2.1-2.3). 
Let us therefore introduce the initial amplitude-and phase noise modulation 
rates m(t) and ¢(t),varying slowly in time and such that m(t) << 1 (small 
amplitude fluctuations) and o(t) << 1/2 (small phase deviations), where 
both m(t) and ¢(t) are assumed to be truncated sample functions of wide- 
sense stationary random processes. Furthermore, let us express the 
perturbing noise voltage near the carrier frequency by means of these 


rates and the output carrier voltage Me as follows: 


v(t) 


Re{V [1 + m(t)]expjlw t + o(t)]}, 


|V.{m(t)cos[wt + arg(V.)] - IV. lo(t)sin[w t + arg(V.)] 3 
(2.20) 
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where the subscript s denotes the stochastic part. Comparing Eqs. (2.20) 
and (2.8) and realizing that Ve has been chosen to be a real number in 


the oscillator model of Fig. 2-la, we obtain 
Vat)emoViemGbeetze av Gt) az Vo(t) : (2521) 


Hence any conclusions regarding the statistical or spectral properties of 
m(t) or $(t) are equally valid for v(t) and v(t)» respectively. 


The Fourier transform of Eq. (2.20) yields 


o) > Jo(wtw.)] 


(2522) 


Velm(w-w,) + J6(w-a,)] + VeLm(wtw 


Since m(t) and o(t) are restricted to low modulation frequencies, the 
upper and lower sidebands of v(t) for a particular modulation frequency 
a << w are 
“gh . 
v(ugta) = 5 Velm(a) + §o(2)] (2.23) 
* * 
v(w_-2) = 5 VeLm'(a) + jo (2)] (2.24) 


The Fourier transform of the initial modulation rates can now be 


determined from Eqs. (2.23) and (2.24) in terms of the sidebands: 


v(wt2)/V, + v(w-2)/V (2.25) 


m(2) 2 


(a) = -§{v(wta)/V, - vi (w-a)/Ve} (2.26) 


(The asterisks in Eqs. (2.24 - 2.26) denote complex conjugate quantities). 
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If the equivalent circuit in Fig. 2-1b is used in place of that in Fig. 
2-1la, v(w) in the above formulas would be replaced by i(w) and Ve by 
Tos i.e., by the fundamental current amplitude through the real part of 
the load admittance. 

Since the quantities V(2) and v, (2) INSEQS. “2.46) sand) (2.14 h.¢an 
be replaced by m(2) and o(2), it also seems proper to normalize v(2) in 
Eq. (2.19) so as to obtain a completely normalized right-hand side of the 
system (2.16), (2.17) and (2.19). This can be done by relating v(2) to 
V4 sthe dc voltage across the active element,Denoting the new quantity by 


by (2) and combining it with Eqs. (2.25) and (2.26), we obtain in matrix 


form 


m(2) V/V. 0 ay v(atw.) 
iC) cn © ey v(2) 
(2) HA EE Swlaaaey | (2.27) 


where the property v(2-w.) = v (u.-2) is used. 


2- cn2ipStat istice Wand ispectral: Properties: of the Initial) Noise 
Modulation Rates and of the Perturbing Noise Voltage under 
Small-and Large Signal Conditions 


Since m(t), o(t) and b(t) are assumed to be truncated sample 
functions of stochastic processes, the power spectral densities Sam (%) » 


S44 (2) and Sip (2) of m(t), #(t) and by(t), respectively, can be calculated 
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from the corresponding Fourier transforms as follows<“°<?: 
S_ (9) = lim E{|m(a)|°/2T3 (2.28) 
mm 
To 
S_ (2) = lim E{]o(a)|°/2T3 (2.29) 
o> Tee 
S.. (2) = lim E{|b.(2)|°/2T} (2.30) 
bb T Vv 
—>0o 


where E denotes the expectation integral and T is the time boundary used 
in the truncation of m(t), #(t) and b(t). Similarly, the cross-spectral 


densities are given by 


. * MES 
ng = a E{m(2)o (2)/2T} = > om (2:03) 
re * meat 
Sib = ae E{m(2)b, (2)/2T} = Spm (2432) 
: * * 
Sb = lim E{$(2)b, (2)/2T} = She (2-933) 


T0 


Substitution of Eqs. (2.27) into Eqs. (2.28) - (2.33) yields expressions 
for the auto-and cross spectra of m(t), o(t) and b(t) in terms of the 


spectra S yyw) oF Vlt)s 


Sam(2) = [Syylweta) + Sw -A) IV, 1? 


+ Gan EtReLv(w,+2)v(w.-2)/V“I/T) (2.34) 


T>00 


iP Pit ; “a 
; Ra 
ry €§ | 


| (88.5) 
(@$.5) 


(0€,$) . ‘ 


beau vats ott off at T bas silat 


fst¥peq2- 22070 odd Inet fake (RG d bra ise tsonn 
oa | on a8 8 ¥ 


(Te. S) 
($€.8) 


(E608) iil 
| poe ee ane 
7 a dy, atl se, in . ia wl ave ‘ah 1 aah noi 


Lac 


‘Oey mt duvet Aa 2 at a : 


ee, “hal 


mt Nal ee -) a7 wie vee 


24 


S.(0) = [5 (uta) +S (w-0)I/1V, 17 


$¢ wie 
- lin E{Re[v(w,+2)v(w.-2)/V.“1/T} (2.35) 
Spp lt) Syl” (2.36) 
Sp (0) = Vin E(Im[v(u +0) v(w,-2)/V-1/T) 
+ I[Syylogt@) = Syylug-@)I/1VgI* (2.37) 
Sp (0) = Lim E{v(w ta)v"(a)/VV. + v (wa) (a)/VVo} (2.38) 
Syp(@) = Lim EC-JLv(upta)v (2)/VoV, - V" (wea) (V/V VB. (2.39) 
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As mentioned in Chapter I, the perturbing voltage v(t) is normally 
postulated to be independent of the amplitude of oscillations. ve is 
then reasonable to assume that v°(t), i.e., v(t) around the frequency Wes 
is a sample function of a wide-sense stationary, normal, narrowband 
noise process. Also, the low and high frequency part of v(t) are assumed 
to be uncorrelated. Alternatively, the total perturbing voltage is often 
assumed to have a white spectrum. All these assumptions imply the 
statistical independence of the Fourier transforms v(w +2), v(w-2) and 
v(2), so that all the cross-spectral terms in Eqs. (2.34) - (2.39) become 
zero. Then 


(a) = [Syylugta) + Syy(wg-2)1/ [Ve |? (2.40) 
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Spl) = Syq(2) = HUSyylugta) - Syylug-@)I/ [Ve (2.41) 


$2 (dna) Stel eau es nu: (2.42) 


Under large-signal conditions, the noise voltage v(t) may be 


amplitude dependent. In view of the existence of nonlinear transformations 


of the basic noise process in the active element, which manifest themselves 


in the form of frequency conversions, the requirement of wide-sense 

stationarity, let alone normality for v(t), appears to be too restrictive. 

Incidentally, even with m(t) and 9(t) wide-sense stationary (these 

assumptions will be kept throughout this study), v°(t) as expressed by 
23 


Eq. (2.20) is generally non-stationary’~. As a result, the two auto 


spectra Sam 62) and S,,(2) are not necessarily equal and the cross spectra 


b> 


Sing (2) and S m6) are not purely imaginary. Furthermore, the spectra 


$ 
Sab (2) and S 5b) do not, in general, vanish. The practical consequence 
of this generalization is that all the auto-and cross spectra must 
generally appear in the final formulas such as those that may be obtained 
from the perturbation equations (2.16), (2.17) and (2.19). This is shown 
in Sec. 2-4. Furthermore, these auto-and cross spectra must be obtained 


by applying the general formulas (2.28) - (2.33), as shown in Sec. 3-4, 


rather than from any simplified relations such as Eqs. (2.40) - (2.42). 


2-3 A Review of the Quasistationary Perturbation Theories of Thaler et al., 


and Viaardingerbroek 
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Thaler et al. ~ oriented their theory towards IMPATT diode oscillators. 


es 


33s) ah. * VAG ade - 


ad yom (3)N Spedfov seton ait « 20} SFB! ghe~s 
enofisimotensys sanilnon Fo Spnaterxe ont to woitw. a _tnabaaaab 9 


eavi sensi t25? tos ni Fe  Sismelp evisoe ett ba ) 
senge-sbTW bef sisi Fup ony. enone hee to 
evivatytean obs od od 21e90Ra) .(t)v IO? VI TiRarOR Sots Mbt: 
seedt) yrenotesde senge~sbew (3)— bas (f)m fttw asvo vette 
vd bee2siqxones Bey ¥ <(ybuae ates suorlouordt Heeaied ffiw 2 
otus ows Sit .tiuesy ap 2A A eysapisete-nn ‘gi fevsnap at (08:9) « 
By339q2 “BFov9, art bas. fevps viiypaegasn ton a6 (eh 2 nw (gg? 6 i “4 
s4J48qe IT Taine vishal sae ish Tommi \isvuG fon ove 40)... ie brs. og 3 
s3neups2naa fsoittonnyy sat suena em (a 
$2uh grtseqe:deoNs a a Mp irs oersipalbadinn is 


bontstdo sd vem! tedg aeod3! aA dave resi nn o% 


Pag i 


nwone ef 2FaT (OTS) bas (nah iat 5) oman ae ot 
: 


penretdo-sd. toun axdange zor j 


‘B-E Jo92 nt sfwonkt 25) «(6E:9) jae Rssilaioal toils oil seat wd i. 
1(Sb-8) 1 (Ob :8) a te es tate mot sas voi 1 


if 


rn 


“yexdtsthis26 sboth TARE! Bangle wont ‘yteddls 3 
ie oe - ae ; | <6 es 


26 


For this reason, they included the effect of bias current fluctuations 
in perturbation equations, which had not been included before in such an 
explicit form. The resulting perturbation equations are of the form of 
Eqs. (2.13) and (2.14) with the frequency dependence of the active 
element impedance neglected. Furthermore, the bias current fluctuations 
AI, (t) were assumed uncorrelated with the perturbing voltage v*(t) near 
the carrier and therefore uncorrelated with the components v(t) and 
v(t). Such an assumption is certainly justifiable at low oscillation 
amplitudes, i.e., in a small-signal theory. As a consequence, no bias 
circuit equation such as Eq. (2.19) is really needed and the system of 
equations (2.13), (2.14) is simply solved for aA‘ (t) and Aw(t) in terms 


of v(t), v,(t) and Al, (t) as the known quantities. In our notation, 


Thaler e¢ al-obtained’® 
dX dR aR. dy axe dR 
C Page drums Eyes aCe tls 
3 vV (ta ‘ yoo ‘ AI (t)A Loy W al W ] 
Mae) =i. C OLNIC OF RAC 
aR” dX) kes dR, 
A\(—+ — - + -—] 
Aue endleecec 
(2.43) 
ae aRe Takata Cole te Oe 
Sie ee eyes ee 
v_(t) ; ] 
Cc 5A S 5 Al b C ar, sl a1, 5 Al 
Aw(t) & C C ; C Cc 
op ORa dX, aX0 dR, 
AL a= - ——] 
. ane ach Mey: mene 


(2.44) 
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To simplify these relations, the following notation is introduced: 


iE ae Ls L : 
=— + ja— = |=—Ilexp(jo) (2.45) 
du dw. dw. 
aR aXe ae 
mee tf joer allie expdin’ (2.46) 
aA. 3A. A 


lexp(jv) . (2.47) 


Furthermore, the small-signal assumptions (2.40) - (2.42) are used so 
that relatively simple expressions result for the noise spectra Spy (Sr) 


and Soy 62) of AA (t) and Aw(t), respectively: 
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wacale 
agi?) S gel) (Ag) FI sin(o-n))-* ay 
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S11(2) (2.48) 


ii ee, ee | 
Ss5tP Ac dpa! S'nien)) + Key S51 (2) ‘ (2.49) 
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where Soc (2) = S., (2) are the power spectral densities of v_(t) and ve (t)> 
respectively, where S112) is the power spectral density of the bias 
current fluctuations AI, (t), and where Kame Key are given by 
yas azn 
Kay = [s72lsin(v-e)/{|—?1sin(e-n)3 (2.50) 
AM ol 1 
0) oA 
c 
aZq dz, 
Key = Fey A LEA eta ! (2.51) 


The second term on the right-hand side of each of Eqs. (2.48) and (2.49) 
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is usually referred to as "up-conversion noise" or "modulation noise", 
the latter being more appropriate for reasons given in Sections 2-4 and 
2-5. Eqs. (2.48) - (2.51) are the fundamental results of the theory of 
Thaler et al. and, in this thesis, we shall not comment on the detailed 
influence of the angles 6, n and y. 


2 introduced a bias circuit equation in a form 


Vlaardingerbroek 
sami tar to Eq. (2.18). e this was done partly in an effort to obtain some 
information on the large-signal behaviour of oscillator noise, since the 
effect of down-conversion by rectification [as represented in Eq. (2.19) 
by the terms proportional to AA‘ (2) and Aw(2)] is likely to take place 
at quite large amplitudes of oscillation. The bias current fluctuations 
AI, (t) are taken as the unknown quantity and the system of Eqs. (2.16), 
(2.17) and (2.19) is solved for the three measurable quantities aA'(2), 


Aw(2) and AT, (2): 


&] is 
aR ae aR” 
aA’ (2) 4 a ot VQ (@)/A, 
aA. C fe) 
ayer dX at, 
Aw (2) = —+ 9 — r) aT -v.(2)/A, FY 
aA. C 0 
a7 Swi 
AI, (2) Lal eo IO) -v(a) (2.52) 
P) GC 


where Z7(2) is the total impedance (diode + external circuit) in the bias 
CIV CUIT. 
The inverse of the matrix in Eq. (2.52) is not given in Ref. 19 for 


the general case. Rather, the inversion is carried out after the intro- 
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duction of the following simplifying assumptions, valid specifically 


for Impatt diodes: 


ned: 360 le $490 gees 

3Al aq a, 3A! ; ‘nie 
(e E 

Q wW W W W 

9Z aR aR® dI aX 9X” dI 

1:0; -—¢-- 53-7 ,-#--i4 (2.54) 
7C aA. =O dA aA. es dat 


Furthermore, the total bias circuit impedance Z>(2) is assumed real, 


Z,(2) = Rr (2). Eq. (2.52) then assumes the form 
[ dx, dR, dR, axe aRe dR, 
ECD aa a ee TD vo(2) 
C Cc qi C ) @) C 
Ry axe aR? 
dw(2) | = =~] - a2 a(ity), spt ality), 0 vlan k 
A_D 0 O 
Cc 
dX dR dP sax” andsR? idx 
L B Epc te hago a ago Folate be 1 
AT. A) Byq—, -BYq—, - B-(q— - =—) A v(2) 
b dw dw. Ry Wo aly aly We 
(2:55) 
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A number of conclusions based on Eq. (2.55) are drawn in Ref. 19. 
One of them is of particular interest to this study; Eq. (2.55) shows that 


the amplitude noise AA! (2) increases as the signal level, and therefore 
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ly|, increases, while the FM noise Aw(2) is independent of the signal 
level if V_(2) and v, (2) are assumed to be amplitude independent. 
According to experiment however“! , the FM noise does depend on the 
amplitude of oscillation; this can be accounted for only if v (2) and 

v, (2) are amplitude dependent and, possibly, cross-correlated. This 
indicates again that a truly large-signal theory must take this RF 
amplitude dependence and cross-correlation into account in some systematic 


manner. 


2-4 _ Large-signal Modification of the Theories of Vlaardingerbroek and 


Thaler 


Vlaardingerbroek |? has not given any explicit expressions for the 
power spectral densities of the quantities in Eq. (2.52). However, his 
approach is implicitly a large-signal one and all that needs to be done 
~ is to calculate the required spectral densities and allow for the RF 
amplitude dependence of, and cross-correlation between, the noise 
voltages Vo(2)s v, (2) and v(@). 

Let us first normalize Eq. (2.52) by substituting into it both the 
output- and initial noise modulation rates. For this reason we will also 


normalize AT, (2) to the dc bias current Ty: 
Br (2) = AT, (2)/1 : (2.57) 


Eq. (2.52) then assumes the form 
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M(2) m(2) 
B(a)} =  [P']  |oy(a)} 
o(2) o(2) (2.58) 


where the matrix P' is the inverse of the perturbation matrix P given 


by 


& 
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Zo 5h! Ly eo 80. ’ 
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Fe nO ne te aes (2.59) 
SES eMac a 
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with Zo = Vol 5 = dc impedance (resistance) of the active element. 
The power spectral densities Sram (2) » Sap (2) and S59) corresponding 
to M(e)., BL (2) and o(2) are given in terms of the initial auto- and 


cross spectra as follows: 
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The coefficients of the inverse matrix p! are, in general, complicated 
expressions and it is not practical to present them here in analytical 
form. If the general formula (2.60) is to be used, it is better to first 
compute the coefficients of P numerically and then to invert the matrix. 

If the effect of ‘enoddlat ion is neglected, as in the study of 
Thaler et ane, considerable simplification can be achieved as far as 
the coefficients of p! are concerned. Yet the final result is still in 
the form of Eq. (2.60) since under large-signal conditions the cross- 
Spectra cannot, in general, be neglected. From the point of view of 
accuracy, it may be expected that a large-signal theory neglecting 
demodulation will be satisfactory up to very large amplitudes where 
demodulation becomes strong. In the case of a large impedance in the bias 
circuit, the demodulation terms in Eq. (2.59) become insignificant and the 
generalized Thaler's theory is then applicable at all amplitudes. 

The effect of demodulation is represented by Poy and P53 in the matrix 
(2.59). If these are set to zero, the matrix can be easily inverted and 
the coefficients for calculating Sym (2) and S 5962) in Eq. (2.60), using the 
definitions (2.45) - (2.47), are 
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The auto- and cross ppeceral densities of the initial modulation 
rates, needed in Eq. (2.60), can be determined from the equations describing 
the active element. This will be shown in Chapter III. Prior to that 
however, a simplified version of the above generalized theory is presented 


in the next section. 


2-5 Simplified Large-signal Version of Thaler's Theory 


The large-signal generalization in the previous section resulted in 
an equation which contains a large number of terms compared with the small- 
Signal approximation. It is almost certain that some of them will be 
relatively unimportant; intuitively one would expect that the amplitude 


dependence of the two auto-spectra Sam) and S, (2) will contribute much 


ox) 
more to the amplitude dependence of Sym (2) and S 59 (2) than will that of the 


cross-spectrum Sin (2). The decisive factor here is the behaviour of the 


> 
active element impedance with respect to RF amplitude and frequency and the 
external circuit properties. 

It turns out that a very simple approximation can be obtained if the 


following restrictions are imposed on the oscillator: 


1) The frequency dependence of the load resistance is negligible compared 
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to that of the load reactance, near the carrier frequency Wo! 


dR, dX) dX) dZ) dX) 
ass dim if so that do ra signa) ; (2.67) 
€ c 


2) The reactance of the active element changes much less with the RF 


amplitude than its negative resistance: 


De ao Ray one Re 
|| <<|—} so that exile - —# : (2.68) 
3A. 3A. 3A. 3A. 


The first restriction is in good agreement with reality if a simple 
external circuit is used. The other restriction is an idealization; 
altogether they characterize the simplest oscillator conceivable. Its 
active element would be represented by an amplitude dependent negative 
resistance with either no reactance at all or with an amplitude 

| independent reactance that can be included in the external circuit. The 
external circuit itself is a series combination of a frequency independent 
resistance and the necessary reactance to form a high Q selective circuit. 


With the above assumptions, the output noise spectra can be written 


as follows: 
ge PL YZ 
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S59 (%) = | pals 26 5 (2) + 2Re(P3P44)S 44 (2) + Paglase Soin (2) ys 2770) 
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is eha eas 
Pl, = {5- | Pl} (2.71) 
11 R 1 
lp vey 
Cc 
dZ 
oy oypaiairy Lara 
3 {da |? (2.72) 


and Pro and oD, are of the same form as in Eqs. (2.63) and (2.64). We 
note that Eqs. (2.69) and (2.70) with the coefficients as given by Eqs. 
(2.63), (2.64), (2.71) and (2.72) are very similar to the small-signal 
equations (2.48) and (2.49). The main difference is that Sam 6) and 
S44 (2) are assumed to be amplitude dependent and are not declared equal. 
Also, the cross-spectral term is missing in Eqs. (2.48) and (2.49). 

The first terms on the right-hand side of Eqs. (2.69) and (2.70) 
may be called "primary noise", whereas the third terms represent modul- 
ation noise. The latter contribution is often called up-conversion 
noise in the small-signal theory. It should be understood however, that 
under large-signal conditions even the primary noise may contain up- 
~ converted components because of the amplitude dependence of the initial 


spectra Sam’) and S,,(2). It is therefore suggested here that the term 


o> 
modulation noise is more suitable for the third terms in Eqs. (2.69) and 
(2 a Oc 


The primary noise contribution in Eq. (2.70) can be also rewritten 


in terms of the quality factor of the external circuit: 


| 2 
Sop (2) = S4,(9)/(Q,, (2/a,)} (2.73) 


The often used "rms frequency deviation" is given by 
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A = (9/2r)[S,,(a)] 7, (2.74) 


rms‘) 
which for primary noise only gives the well-known formula 


2 


afrms (2) = (F./0,, JES, ,(2)] (2.75) 


rms 

In summary, in the simple oscillator characterized by the conditions 
(2.67) and (2.68), the internal amplitude-to-phase or phase~to-amplitude 
conversion, as indicated by the diagonal lines in Fig. 2-2, is so small 
that the primary amplitude (phase) noise at the output can be expressed 
in terms of the initial amplitude (phase) noise only. 

The above approach yields the simplest large-signal approximation of 
oscillator noise. The reason for introducing and discussing this simple 
model rests in the fact that the general model, described by Eq. (2.60), 
may sometimes reduce to the former even if the condition (2.68) is not 
satisfied. This depends on the relative magnitudes of the initial spectra. 
This fact is utilized in some of the calculations of the output noise 


spectra presented in Chapter V. 
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CHAPTER ITI 


SPECTRA OF INITIAL (PRE-OSCILLATION) NOISE MODULATION RATES: LARGE SIGNAL 
CALCULATION 


In the previous chapter, formulas have been derived that employ 
the spectral densities of initial noise modulation rates or, alternatively, 
the spectral densities of the two orthogonal noise voltage components 
v(t) and v,(t) in the representation of the perturbing noise voltage. 
As shown in Sec. 2-2-1, these quantities can be determined from the side- 
bands of the perturbing noise voltage v(t). A general procedure for 
computing the latter quantities and the initial modulation rates will now 


be given. 
3-1 Excitation of the Active Element 


The equivalent oscillator circuits in Figs. 2-la and 2-1b suggest 
that the noise source voltage or current could be computed from the 
appropriate equations of the active element which is either open- or 
short circuited. Under small-signal conditions, the noise mechanism can 
be considered "classical", i.€., independent of the RF voltages or 
currents. in the circuit, and only the dc level needs to be considered. 
Under large-signal conditions, we must find a way of simulating the 
effect of large RF currents and voltages present at the active element 
when the circuit is oscillating. This can be achieved by driving the 


active element from a deterministic current source (if the open-circuit 
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noise voltage v(t) is to be determined), or by a voltage source [if 

the short-circuit noise current i(t) is to be determined] at that 

amplitude and frequency at which it would normally oscillate. Circuits 
corresponding to this desired situation are shown in Figs. 3-la and 3-1b. 
Note that they are non-oscillatory circuits. In singly-tuned oscillators 
with quality factors that are not too small, either the current or the 
voltage at higher harmonics of w, will be relatively small and one of the 
circuits in Fig. 3-1 will be applicabte. Clearly, the circuit in Fig. 

3-1a (3-1b) must be used in conjunction with the oscillator model in Fig. 
2-la (2-1b). The subscripts s and d refer to stochastic and deterministic 
quantities, respectively. 

One point worth mentioning is that the bias circuits in Figs. 3-la 

and 3-1b are not open- or short circuited. This is so because we wish to 
take into account the possible effect of bias circuit impedance on the 

open-circuit noise voltage (or short-circuit noise current) at high 
frequencies. This effect, if it exists, might be called "initial modulation" 
as distinct from the modulation mechanism mentioned in the previous chapter. 


As a side note, it is suggested that the often encountered term "up- 
conversion" should perhaps be used in a broader sense than is typically 
accepted; namely, to include modulation, initial modulation and the conversion 
which takes place even if the bias circuit in Fig. 3-la (3-1b) is open 
(short)- circuited. This latter effect is due to the mutual correlation 


between the open-circuit noise voltages (short-circuit noise currents) 


at various frequencies, under large-signal excitation. 
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Ig=1) cos We t 


V(t) = V(t) + Vd 


fe 8) 
Vg 21/2 » Vnexp (jn We t) 
-@ 
I(t)=I}cos wet + Ip(t) 
Fig. 3-la EXCITATION OF THE ACTIVE ELEMENT BY CURRENT SOURCE 
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V(t) = V) cos Wet + Vp (t) 


Fig. 3-lb EXCITATION OF THE ACTIVE ELEMENT BY VOLTAGE SOURCE 
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3-2 Mathematical Description of the Active Element; Generalized 
Langevin Equation for Stochastic Quantities 


Let us assume that the nonlinear element can be described by the 


differential equation 
BT 
ay{1(t) V(t) 3 + bgCt(e) WO] SH + a, c1(e), wid SAD + 
i 


+... tay[1(t).v(t)] SE) + by or(t).v(ty] ME) + c orcty.w(ty] = F(t) 


(371) 
where either I(t) or V(t) is a given periodic driving signal as discussed 
previously and F(t) represents the internal noise mechanism in the active 
element (thermal noise, random ionization,etc. ). The statistical and 
spectral properties of F(t) may, in general, be dependent on some 
controlling parameter related to the power level, i.e., either to the 
% ie 


average of |V(t)|~ or of |I(t) 


Under non-oscillatory conditions, it is reasonable to write 


| | 
—_— 
ct 
~— 
iT] 


I ,(t) + I. (t) T(t) << I(t) (ena) 
V(t) = V y(t) ul V(t) Vv (t) << V g(t) 3 (3.3) 
where V y(t) and V(t) (1,(t) and I.(t)) are the deterministic and 


stochastic components of voltage (current), respectively. It may sometimes 


be convenient or even necessary to formulate Eq. (3.1) in terms of other 
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electric quantities (e.g. electric field) or to include the latter in 

Eq. (3.1) in addition to voltage and current. This does not basically 
change the procedure outlined here; the linearization in stochastic 
quantities indicated in Eqs. (3.2), (3.3) is simply carried out for all 
time dependent variables in Eq. (3.1). Eq. (3.1) can then be decomposed 
into two equations; 

1) deterministic, which is of the same form as Eq. (3.1) except that the 
right-hand side is zero 


2) stochastic, with the second- and higher order products neglected: 


drt) 28, 24, dT 4(t) 
a LI 4 ( ) V4(t)] Seagrass Riese te UL aaa ne b Li g(t) V4(t)] 
dt d ate 
a Welty. 3b ab dy. (t) aC 
S m m d 0 
X + [—— I(t) + = V(t)] + + —— I(t) +t 
atm ory S avy S at™ oly S 
8C, 
* ayy Volt) = F(t) fea 


This iS a generalized version of the Langevin equation. 


3-3 Fourier Transform of Langevin Equation; Open-Circuit Noise Voltage 
(Short-Circuit Noise Current) 


The Langevin equation (3.4) is, in general, a nonlinear differential 


equation since the coefficients an? band their derivatives with respect 


m 
to I, or Va are functions of the latter and therefore general functions 
of time. However, the specific excitation of the active element, as 


introduced in Section 3-1, enables us to simplify the above-mentioned 
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equation to a linear differential equation with periodically time-varying 
coefficients. Since one of the quantities Ty and Va is enforced as a 
periodic, sinusoidal quantity, the other may be approximated by Fourier 
series aS indicated in Figs. 3-la and 3-lb. Thus we have 


either 


tH 
i} 


ten 
q = 1). cosu.t 3 Mi = gen exp(jnw t) : (3.5) 


or 


=< 
iT] 


baal fe 
ae cosw.t 3 i 3 ) I exp(jnwt) : (3.6) 


The coefficients an? ba and their derivatives are then also 


approximated by Fourier series 
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a: = bam nexP(dnwct) pe Dra= seb nex? (ina, t) ; (337) 
ete. 
The coefficients ares b, no? cre are functions of the harmonic components 


of voltage V4 and current Ty: 
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or 


men br mV 22] ty slosd_ps-- +) ; (3.11) 


etc. 
The substitution of Eqs. (3.5) - (3.11) into Eq.. (3.4) leads to an 


equation of the form 


) ee : ; erst) 
C. exp(jnw _t )————- + D. exp(jnw_t )—————_ 
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(35.12)) 


On making the assumption that I(t), V(t) and F(t) are truncated 
sample functions of a random process, the Fourier transform can be 


applied, resulting in an equation of the form 


‘ AV, (wtnw) a BL, (wtnw.)= F(w) : (3 ais) 


N=-o 


where the coefficients Ans B. are functions of w, win and of the Fourier 


n 


components was Vy Cra I, brad Sa of the deterministic voltage and 


Ne) Os 


current. The latter can be determined from a separate large-signal 
analysis using Eq. (3.1), where F(t) = 0, in conjunction with the 
excitation diagrams in Figs. 3-1. V(wtnw), I (wtnw ) and F(w) are the 


Fourier transforms of V(t); I(t) and F(t), respectively, at the 
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frequencies indicated. 


Since we are interested only in the spectra at the frequencies 


w=Q+ &w |al<< ws 2 = integer , (3. 
C @ 


we can substitute Eq. (3.14) into Eq. (3.13), and obtain 


5 A (atew )V [at (nt2)u J+B, (atew )I [ot (nt2)w J =F(2+ ww), (3. 


N=-o 
which,introducing p = n+, can be rewritten in the more useful form 


M 


aes 


2,2 
p= P : 


because in actual computations, p will be limited to a finite range. 


Incorporating the relations (c.f. Figs. 3-1) 


Vg(2) = Z(a)I,(a) , Mite: 
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we obtain a system of linear equations of the form 
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or a similar system for current. An inversion of Eq. (3.19) yields, in 
matrix form 


j 


[V.p] = [cy 


Je glithyl : VE5 = V.(atpw.), Piri F(otLu, ) va peli3ia20) 
If the above procedure is oriented so as to obtain Eqs. (3.19) and (3.20) 
in terms of a stochastic quantity other than voltage or current, the 
corresponding relation between the quantity used and the open circuit 
voltage (or short-circuit current) must be substituted into Eq. (3.20) to 
obtain the latter. This will only modify the coefficients Crnge a Then, 
for a particular p = + q we obtain the upper- and lower sidebands 
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g=-M 
ae 
V,(2-qu..) = ai C gee gah (at 2H.) ; (3.22) 


and for the noise voltage at baseband frequencies 
Mites 
VA(Q) =a04 F(at£u) 4 (3.23) 


The form of Eqs. (3.21), (3.22) and (3.23) clearly indicates potential 


conversion in the active element. 
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3-4 Fourier Transform and Power Spectral Densities of the Initial 


Modulation Rates 

It will now be postulated that the Fourier transform of the perturbing 
noise voltage v’(t) in the oscillator model in Fig. 2-la is approximately 
equal to that of V(t) in the vicinity of We 

v(w +2) = \ (w +2) 5 (S524) 
Similarly, the Fourier transform of the perturbing noise voltage v(t) 
is now postulated to be approximately equal to that of V(t) at baseband 
frequencies: 


Wie (Oe. (3925) 


Hence, the Fourier transforms of the initial noise modulation rates at 


frequency 2, from Eqs. (2.27) and (3.21), (3.22), are 


Mo 
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4(2) Cer : Cy g/Ve)F(atan,) (3.28) 


In most cases,F(t) in Eq. (3.1) will belong to a wide-sense stationary 


process with an almost white spectrum so that, on applying Eqs. (2.28) - 
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(2.30) 5 


Lim E{F (atmo )F (atnw.)/2T} rumen ne oe (3.29) 


T30 


The power spectral densities of m(t), by (t) and $(t) are then 


M ] 
7 1 Z 
Sam (2) = gle eategl = + sod 2 JN. | Se p(2tew) (3.30) 
eer 2 
Sip (2) = ey on a8! “0! Spp(2t2w) (3.31) 
ee 
S44 (2) = ey o1-2 0 Me 7 Bae Q, lV | as fF (2+ kw) . (342) 


where Sep (2) is the power spectral density of F(t). 


In a similar way, the cross spectra S_ (2), Sp (2) and Sp () can be 


mo 
determined by applying Eqs. (2.31) - (2.33) on Eqs. (3.26) - (3.28), 


using the property (3.29): 
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Finally, it may be useful to compute the spectra of the open-circuit 
noise voltage at the sidebands. From Eqs. (3.21), (3.22), (3.24) and 
(3.29) we have 


ona! 
if i YA 
Syy (mo#2) = ey 1-240! Spe (2+ hu) (3.36) 
Seacenpeebuie ; esha scan ona (3.37) 
wc ee ce get eel 


where Syy(w) denotes the power spectral density of the open-circuit noise 
voltage v(t) at the frequencies indicated. 

It is seen that all the auto- and cross spectral densities required 
in the general expression (2.60) for the amplitude and phase spectra of 
an oscillator, can be obtained by means of the power spectrum of the basic 
noise process and the coefficients in Eq. (3.20), which characterizes 


the frequency response of the active element to its basic noise process. 


3-5. Computation of Oscillator Noise: General Procedure 


The derivations in the previous sections and in Chapter II resulted 
in a considerable number of equations as well as of expressions to be 
substituted into those equations, to obtain the noise at the oscillator 
output.’ The purpose of this short section is to consolidate the 
fundamental results and to show the procedure for computing the oscillator 
output spectra, given the spectrum of the basic noise process in the 
active element. 
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“interaction diagram" of Fig. 2-2, since it is now known how the perturbing 
noise voltage is obtained from the basic noise process in the active 
element. Such a complete conversion diagram is shown in Fig. 3-2. The 
basic noise process first passes through the stage of initial conversion 
caused by the nonlinearity of the active element and the presence of an 
RF signal. The resulting noise voltage and the corresponding initial 
noise modulation rates are then again subject to conversion when the active 
device is allowed to oscillate freely. 

It is immediately seen from Fig. 3-2 that the simplest and most 
elegant method for calculating the output spectra is as follows. Firstly 
the Fourier transforms of the output spectra are obtained from the Fourier 


transform of the basic noise process via all the matrices involved: 
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The output spectra are then calculated from the following formulas: 
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Sails lin E{M(a)M (2)/2T} = J IR gl Spp(o+2e) (3.39) 
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S59 (2) = Ly E{o(2)® (2)/2T} = Pais rt Spe (ate) ; (3.40) 


where the property (3.29) is used. 

It may sometimes be desirable to know the initial spectra along 
with the output spectra. In this case, Eqs. (3.30) - (3.35) are used to 
compute all the initial auto- and cross spectra, which in turn can then 


be inserted into Eq. (2.60) to obtain the output spectra. 
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CHAPTER IV 
APPLICATION OF THE NOISE THEORY TO IMPATT DIODE OSCILLATORS 


The large-signal theory of oscillator noise, developed in the 
previous chapters in a general form, will now be applied to IMPATT diode 
oscillators. Two decisions must be made at this stage in order to be 
able to achieve a reasonable compromise between simplicity of application 
and sufficient accuracy of the results. First, a suitable mathematical 
model of the IMPATT diode must be selected so that the procedure for 
calculating the initial spectra (c.f. Chapter III) can be carried out. 

The initial spectra are of particular interest in this study because of 
the possibility of comparing their large-signal behaviour with the usual 
small-signal assumptions imposed on the perturbing noise source. Having 
obtained the initial spectra, the modified quasistationary perturbation 
theory described in Chapter II is used for the computation of the measurable 
output noise spectra. Depending on the accuracy desired and on the 
associated complexity (i.e., the tolerable computer expense), this may be 
obtained by using either the most general perturbation matrix in Sec. 2-4 
or one of the two simpler approximations in Sections 2-4 and 2-5. The 
problem of computing time is very real here because of the number of time- 
consuming auxiliary calculations and the large number of combinations of 


the parameters involved in computing both the initial and output spectra. 
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4-1 IMPATT Diode Model 


4-]-] Avalanchin -n Junction 


In the reverse biased p-n junction diode, very small reverse 
saturation currents of electrons and holes flow across the junction. As 
the reverse voltage is increased, the maximum kinetic energy of the 
particles carrying this saturation current increases also. When a 
certain critical electric field intensity has been reached near the junction, 
the electrons and holes become capable of producing other electron-hole 
pairs by impact ionization and the phenomenon of rapid carrier multiplication, 
called avalanche multiplication and/or breakdown ensues. This effect can 
be controlled and even utilized for high-frequency generation. The avalanche 
buildup has associated with it a certain delay time; further, after the 
electrons and holes have passed through the region near the junction, an 
. additional transit delay time is introduced as the carriers drift through 
the remaining part of the diode. Thus, negative resistance can be produced 
at certain frequencies, typically in the microwave region. The diodes 
utilizing this phenomenon are called IMPATT diodes, an acronym IMPact 
ionization Avalanche Transit Time. 

In general, the fundamental equations describing the avalanche process 
form a rather complex system and the existing analytical solutions have 
been obtained by making various simplifying assumptions. From these 
approaches the most general deterministic (i.e., noise excluding) solution 
known to this author is that of Kuvas and Lee". The differential equation 


which they obtained can be considered to be a generalization of the results 
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7a) 26 


of Read-~ and Lee et azZ., ~ which were developed for a device somewhat 
more complicated in structure, but easier to handle mathematically than 
a general p-n junction. This structure, the Read diode, has been used 
almost exclusively in various small- and large signal analyses of IMPATT 
diodes, as well as in noise studies of the latter. Although there exist 
various structures, the majority of devices used in actual practice are 
simple one-sided abrupt junctions, for which the Read diode structure is 
the closest approximation that lends itself to simple mathematical 


description. We shall, therefore, choose it as the mathematical model of 


the IMPATT diode. 


4-1-2 The Read Diode 


- is shown in Fig. 4-1 with 


This structure, first suggested by Read 
its doping profile and electric field distribution. Its most important 
feature is the division of the diode into two regions. In the high-field 
region, called the avalanche region, intensive avalanche multiplication takes 
place when a certain electric field is reached, whereas in the lower-field 
region, called the drift region, the injected carriers drift across under 
the influence of the electric field. 

The fundamental quantities for the mathematical description of 
avalanche breakdown are the so called electronand hole ionization rates 
Ot, and a and the electron and hole multiplication factors My and Mo» 
respectively. The former are defined as the number of electron-hole pairs 


generated by an electron (hole) per unit distance. They are strongly 


dependent on the electric field and this dependence is often approximated 


by 
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where Ot » ES are constants and m = 3-9, or by 


a, or a, = Aexp[-(b/E)] (4.2) 
with the values of A, b and m given in tables for various semiconductor 
materials (typically Germanium, Silicon and Gallium arsenide). 

The multiplication factors of electrons and holes are defined as 


follows: 


He EL (OL Cremer AON TAN)” (4.3) 
where 1 (Ww) and 1 (w) are the hole and electron currents, respectively, at 
the boundaries of the depletion layer of total width w. 

The differential equation for the generated avalanche current has 


been derived first by Read@”. 


The most important assumptions were equal 

ionization rates and drift velocities for electrons and holes and neglect 
of space charge in the avalanche region. Lee et see generalized Read's 
analysis to different ionization rates and drift velocities while keeping 


the same form of differential equation; the latter can be written as 


follows 
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where I(t) is the conduction (particle) current at x = We i.e., the 
current flowing out of the avalanche region and consisting of the 


(t) 


electron current I, ft) and of the hole current lop 


i ae: (4.5) 


and where I. is the total reverse saturation current consisting of the 


Saturation current of electrons, am and holes, I.) 


The quantities M and t in Eq. 4.4 are the "equivalent" multiplication 


factor and the so called intrinsic response time, respectively, which are 


given by the following relations<?: 


UE neiel agat gil. ¢ (4.7) 
1./t = A es - Ioy/Tp (4.8) 
W W 
. a a 
M, = ] mk open (- (a -a_)dx')dx (4.9) 
0 X 
W W 
= cn (-} (a -a_ )dx') ax (4.10) 
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In the above expressions, a (v,) is the electron (hole) velocity and 
ce (c)) is the electron (hole) intrinsic response time. 

In addition to Eq. 4.4, another relation is needed to account for 
the influence of the displacement current through the device: 


dE AC) 
Go) ead Ro Gone aaa 9 ae ‘ (4.12) 


where Ty is the external (terminal) current, « is the dielectric 
permittivity, A is the cross-sectional area of the diode and E,is the 
electric field in the avalanche region. 

In the noise problem, it has been generally assumed in the 
literature that the strongest contribution to the overall diode noise 
comes from the random ionization in the avalanche region. Other sources 
such aS generation-recombination noise, flicker noise and noise due to 
the imperfections in contacts or the junction itself have been neglected. 
The first noise analysis of the dc-excited diode was published by Tager®, 
in which a driving noise source term was added to the right-hand side of 
Eq. (4.4) and its spectral properties were derived from the statistical 


9510511 followed, 


considerations of the avalanche process. Other analyses 
with some modifications, but all under the condition of the pure dc excit- 
ation of the diode. A more general derivation was given recently by 
Kuvas@?: The resulting equation for the conduction current is of the form 
of Eq. (4.4) with an added stochastic term whose spectral properties were 


derived for a combined dc and ac excitation of the diode: 


[a/at + (Mc)7"]I,(t) 2(t, + Flt}ir | (4.13) 
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The statistical average of F(t) is given by@? 

<F(t)F(t+o)> ed(o)I (t) a Qe G2 =G 

errr ora LAN a Ves a Gy ker LD ee Ls (4.14) 
2 2 a 2a_M 2 
T iad n pon 20M, 


Here e is the electronic charge and 6(c) is the Dirac delta function. As 
can be seen from Eqs. (4.9) and (4.1) or (4.2), M, is a function of electric 
field. Under the combined dc and ac excitation of the diode, M, can be 
expanded in Fourier series“?. Also, I(t) as a large-signal solution of 

the deterministic equation (4.4) (see Appendix-B) can be inserted into 

Eq. (4.14) and, after taking the time average of Eq. (4.14) and Fourier 


transforming the result, the power spectral density Se (w) of F(t) has 


the form? 
SL Atay inedt JF ey Swe tata a, = ; , 
ere a (outs Ds) o(ulaK de 4 
2 2 TT 4 Gan + 2 
T 1E np 


where I) is the dc (bias) current through the diode, We is the fundamental 
frequency of oscillation and D1 »K, come from the expansions of My described 
in Ref. 20. Since D, is proportional to Ey; i.e., to the electric field 
component at the fundamental frequency, the second term in parenthesis on 
the right-hand side of Eq. (4.15) represents a correction due to the 
presence of an ac signal. 

According to Kuvas, this correction is small for reasonable signal 
levels so that 
‘ (4.16) 
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which for hast Eoaiod (iGo Me aM a, = Oy» etc.) reduces to Tager's 


result® 


Sep(w) * el ats. (4.17) 


Eqs. (4.13), (4.16) or (4.17) and (4.12) will be used for the 


application of the analysis presented in Chapter III. 


4-1-3 Approximation of Real IMPATT Diodes by Read Diode Equations 


As mentioned in Sec. 4-1-1, the majority of actual IMPATT devices 
are one-sided p-n junctions and the description by Eqs. (4.12) and (4.13) 
is only a coarse approximation for these. Yet many reasonable results 
have been obtained in the computation of various IMPATT diode properties 
with even the simplest form of Eq. (4.13), namely that resulting from the 
assumptions of equal ionization rates and drift velocities for both 
electrons and holes. Employing these assumptions in Eqs. (4.9) - (4.11) 
and substituting into Eq. (4.13), we obtain 


7. OW bt) we 
eet (ef ol ER) aRtton )sebrgre! F(t )S tg (4.18) 
0 


where Ts is the transit time across the avalanche region: 


sae wo/V ‘ (4.19) 


Eq. (4.18) is the original Read equation and although iu. 1s strictly valid 


-only for Germanium, it has been used for the analysis of Silicon diodes, 
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which are most commonly used in practice and for which the ionization 
rates Ot and Oy are different. In the latter case, the ionization 


integral in Eq. (4.18) can be approximated as follows?: 


Ww W 

a 4.5 
if Ae | aay (4.20) 
(@) (6) : 


a = (a,-o,)/In(a,/a,) * (4621) 


and where the average ionization coefficient a is assumed to have a 
similar dependence on the electric field as Oy or a, TPE dSe (4% 1) or 


(4.2). Specifically, we shall later employ the approximation 


= Beart] eee cl Gi ' 
a(E.) 7 Wa (E./E,) Ws [i+m(E. Baeen P (azo) 
where EY is the electric field in the avalanche region, at breakdown. 

Employing the usual assumption of ES being independent of position in the 


narrow avalanche region, the ionization integral is simply 
ye 
f a(E,)dx = a(E,)w, + 1+m(E,-E,)/E, (4.23) 
0 


As pointed out in Ref. 24, the intrinsic response time +t generally 
differs from the value 1/2, which in a realistic diode must be corrected 
to account for the influence of space charge in the avalanche region and 
that of the drift region. These corrections depend on the ratio of 
ayant coefficients and drift velocities for electrons and holes ina 


complicated manner“. The corresponding formulas would not be practical to 
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use for our purpose; instead, Eq. (4.18) will be rewritten as follows 
(also using Eq. (4.23)): 


dI (t) 


errr yd SIM Eustis a chet eal tall (4.24) 


kt 
a (e 


and the correction factor k will be calculated from the data given by Kuvas9, 
which are based on measurements of representative diodes at 6 GHz. These 
data are summarized in Table 4-1, Sec. 4-3-1. 

Also, if Eq. (4.24) is to be used for a p-n junction, the effective 
avalanche region width must be evaluated from the electric field distri- 
bution?’ #28, Table 4-1 gives Kuvas' values for Wa and also w, the total 
diode length. 

Finally, for one-sided abrupt p-n junctions, the following relation 
holds between the voltage Vb across the diode at breakdown and the maximum 


~ electric field in the avalanche region Bie ps 


are lw foe a) (4.25) 


Eqs. (4.12), (4.17), (4.24) together with the data in Table 4-1 and Eq. 
(4.25) are believed to represent typical IMPATT diodes and will be used in 


the actual computations. 


4-2 Derivation of the Open-Circuit Noise Voltage 


In order to properly apply the mathematical procedure described in 
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Chapters II and III, the appropriate equivalent circuit must be chosen 
for the specific oscillator under investigation. Experiments have shown 
that, in singly-tuned IMPATT diode oscillators, negligible power is 
delivered to the external circuit at higher harmonics, and Kuvas@? gave 

a simple theoretical proof that the external current components at higher 
harmonics of the oscillation frequency can be neglected. This implies 
that the equivalent circuits of Fig. 2-la and 3-la are suitable models; 
the mathematical procedure will be oriented towards obtaining an equation 
of the form (3.20), for the open-circuit noise voltage rather than 
short-circuit noise current. 

Furthermore, the deterministic large-signal analysis which gives the 
dc and periodic voltages and currents and is necessary in the actual 
evaluation of the coefficients in Eq. (3.20) is also affected by this 
choice of the equivalent circuit. The external deterministic current is 
assumed to consist only of a dc component plus a periodic one at frequency 
w_, the frequency of oscillation. The results of the corresponding large- 
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Signal theory of IMPATT diodes, as derived by Kuvas -, are given in 


Appendix-B. 


4-2-2 Langevin Equation for Stochastic Quantities 


If Eqs. (4.12) and (4.22) are substituted into Eq. (4.24), a non- 


linear differential equation results, similar in form to Eq. (3.1): 
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Reece E(t) (4.27) 


(ete elapse l(t) << To a(t) (4.28) 


exd exs exd 


as in Eqs. (3.2) and (3.3), the stochastic part of Eq. (4.26) becomes 


2 
dB (2) dt. (4) E.(t) dE g(t) 
kt. [-eA 2 + = m = [ea - Dedcb 
E..(t)-E dE". ('t) 
+m rye [eA 85-1, (t)] = F(t), (4.29) 


where the stochastic part of the reverse saturation current I. is neglected, 


Since it is small compared to F(t)°. 


4-2-3 Fourier Transform of Langevin Equation; Open-Circuit Noise 
Voltage 


The deterministic quantities La and Ed in Eq. (4.29) will now be 


expanded into Fourier series as follows: 


Tayq(t) r I, a re => cs Lay (dnwt) ; (4.30) 
where 

Lage 1/0) en lites, | (Aeon) 
and 

Eaq(t) = z } E.nexp(jnwt) (4.32) 


(tS,4), Wha > oes .* i hus the: 
im a. or 
(A0)- , ‘leas ec Ma (linc ry wih yh | 
| yale rahi aia 


| sr 
etioned (88.8). p2. 40 169 tee able pate 


oh) | a a : -_ 
i yy 7 


pelo, Ttuoy = 1eq0 | ene 


” 


sd won Titw(@S.8) «pam 


(O€. 1) 


(fe.») 7) : s : | . "Tenino 0 ys! : : on 7 


65 


Similar expansions can be written for the conduction current I g(t) and 


for the voltage V(t): 
) I exp (jnwt) (4.33) 
4 V exp (inw t) 2 (4.34) 


The substitution of Eqs. (4.30) and (4.32) into Eq. (4.29) leads 


to an equation similar in form to Eq. (3.12): 


dE _(t) dE. g(t) 
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Applying Fourier transform to Eq. (4.35) and limiting the number of terms 


yields an equation of the form 
A_E__(wtnw.) + BI (wtnw.) = F(w) , (4.36) 


which corresponds to Eq. (3.13) in Chapter III. 


The coefficients A, and Bae from Eq. (4.35), are 


Ao = JweALm(E, .-E,)/E,-Jwkt, J-mI /E, (4537) 
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Bo = jukt, - m(E 507 EL) /Ey (4.38) 
J * 
ag = (jumeA/2E, )E. - (m/2E, )I. (4.39) 
a * 

BL = -(m/2E JE, - (4.40) 
n#o 

Substituting 

w = Othues | Q| << py Ris integer (4.41) 


(as in Sec. 3-3) into Eq. (4.36) and introducing p = 2tn, we obtain 


M 

pe-M Ap-a gE as (2tPHe) + Bain pleys(atPw.) * F(at2w,) (4.42) 
for, -M < 2 < M. 
Here 

Ano 42 = Anne (w=arew.) 3 Bo H2 42 = Boag (u=at2u,) “ (4.43) 


Eq. (4.42) is the equivalent of Eq. (3.16) in Sec. 3-3. 

The field and current components oa and tone respectively, in 
Eqs. (4.37) - (4.40) must be determined from a separate large-signal 
analysis. In Appendix-B, the Fourier components of the conduction 
current I g(t) can be calculated from the normalized electric field 
component E/E, at the fundamental frequency Wee It therefore appears 


advantageous to replace E- and yon in Eqs. (4.37) - (4.40) by the 
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conduction current components ia The appropriate relations can be 
derived from Eq. (4.12). This equation also enables us to express 
Toys w) in terms of Es () so that Eq. (4.42) may be written in the 
form 


M 


p=-M Ayn gE ag (@tPH.)/E, + F (tw) /mI , 


M<g<M. (4.44) 
The above-mentioned relations are given in Appendix-A. The final 
normalized form of the coefficients Qn-0,2 is derived in Appendix-C. 
Since we are interested in the open-circuit noise voltage rather 
than the electric noise field, the corresponding relation between these 
two quantities must be known; this is also derived in Appendix-A. This 


relation may be inserted into the matrix system (4.44) to obtain a 


similar system with modified coefficients 


Vo(atpw)/V, 2 F(at2w.) /ml 
-M<2 <M (4.45) 


from which V,(atpw)/Vp is obtained by inversion. An alternative is to 
invert Eq. (4.44) first and then. compute the open-circuit noise voltage 
from the electric noise field. These two alternatives may not be 


equivalent from a computational point of view, as will be mentioned in 
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the next section. 
Once the relation for the Fourier transform of the open-circuit 
noise voltage is known, the procedure described in Sections 3-4 and 3-5 


is applied to compute the initial and output noise spectra. 


4-3 Computed Noise Spectra 


4-3-1 Computation: Conditions, Input Data, Procedure 


The formulas derived in the previous sections are used to compute 
both the 'initial' (pre-oscillation) and output modulation spectra of 
IMPATT oscillators employing a simple series resonant circuit. The 
computations are carried out for Si, Ge, and GaAs abrupt junction diodes 
with physical parameters as given in Table 4-1 for various values of the 
electric parameters specifying the operating point. The physical para- 
meters themselves were also varied. This is important, since they are 
only equivalent parameters that originate, for example, from an effort to 
fit the Read diode equations to an abrupt junction diode. As such, they 
are subject to considerable variation as evidenced by the differing data 


8,10,20,28 Besides, actual diodes 


deemed typical found in the literature 
do differ in their structures even if designed for the same frequency. 
In the calculations the series resistance of the diode is assumed to be 
included in the external circuit and its effect is taken into account 
only in the final comparisons of computed and measured spectra, based on 


the output power rather than diode voltage. No noise sources other than 


that associated with the random ionization in the avalanche region are 
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considered. 

Since the initial spectra are of considerable interest in this case, 
the "all-matrix" method suggested in Sec. 3-5 is not used; instead, the 
initial spectral densities are calculated from the corresponding 
formulas as given in Sec. 3-4 and then the computation of the output 
Spectra is carried out. In this way, a better insight is obtained into 
the mechanism of noise increase in IMPATT diodes under large-signal 
conditions. 

We shall now describe the computational procedure step by step, 
bringing together all the key equations in the proper order and hopefully 
making the application of the method described in this thesis clear. The 
procedure makes use of the normalized quantities introduced in Appendices 
-A and B. 

Calculation of the initial spectra: 

Step 1: The following parameters (introduced in Appendix-A) characterizing 


an IMPATT diode, are assigned specific values: 


2 
m,(w./w) ok wT qsWa/WqsWy/Wo- 


From these m and k are fixed as found in Table 4-1, and WoT and w,/Wq are 
varied around the values given in Table 4-1. The parameter (w/a)? is 


varied in the range 0.1 - 0.7, which is typical for IMPATT diode oscillators. 


The parameter (w/o) is typically in the range 10; * Gud 0 Fe 
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TABLE 4-1 SOME PARAMETERS OF REPRESENTATIVE IMPATT DIODES 


AT 6 GHz (From Ref. 20) 


Step 2: For a given combination of the above parameters, several 
consecutive values of the normalized electric field component 
le.7| = JE 7 /E,| at the fundamental frequency w are chosen, typically 
in the range 0.02 - 0.3. This covers the range of electric field 
magnitudes occuring in practice. To each of these there corresponds a 
certain conduction current (normalized with respect to the dc current), 
whose components are computed from (see Appendix-B, Eqs. (B.2) - (B.5)) 
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where 


i106. 4)/12 (Bess) (4.47) 


with I. and ie being the modified Bessel functions, and 
a -] 
D = (m/k)(woty w./Wq) ; : (4.48) 


The number of components needed must be determined by trial and error, 
until the result in Step 6 does not change appreciably with the increase 
in the number of components. It turns out that 5 harmonic components 
are sufficient for this purpose. 

Step 3: The normalized diode voltage component v1 at frequency We is 


computed from (see Appendix-B, Eq. (B.8)) 


Vy ig felw)-1] 
Vi thayat tl ee =J) e018 ‘ (4.49) 
b Jj(w,/w.) “D(1+w. /wq) 
where 
6(w) = [1-exp(-ju.t4)]/iwoty ; (4.50) 


The normalized output voltage component vA at frequency Wo is then computed 


from (see Appendix-B, Eq. (B.10)) 
Mace V/V, = Re (v1) . (4.51) 


The normalized external current component Pat at frequency We is determined 
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from (see Appendix-B, Eq. (B.7)) 
ie ey T at eu Dy Mo) cle 1 (4.52) 
ex] exl’ “o cl ac al ? ; 


and the normalized impedance Z(w/) of the external circuit, needed to 


sustain oscillations, is given by (see Appendix-B, Eq. (B.9)) 


z(w,) = w CyZ(w,) = (D/2) (w/a) Vy liggg (4.53) 


Step 4: The normalized baseband frequency parameter Of, is assigned a 


8 4 


value, typically in the range 10 ~ to 10 ° and the normalized impedance 


of the external circuit at baseband frequencies, 
z(a2) = wC4Z(2) A (4.54) 


is assigned a value in the range 0 to »~. 


The coefficients An -0 .2 in the conversion matrix (4.44), 
i E__(a2tpw_) M F(Q+2w_) 
te Ee e__(atpw.) = 
Poem pes o£ E. oy, Pre 2k as © mI, : 


(4.55) 


are calculated from Appendix-C, Eqs. (C.13), (C.14), (C€.15) and (C.16): 


For p-2 = 0 


(2 oy(2) E+ a2 +2) IEv(arew,)-1] . (4.58) 
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For p-2 # 0 


1-y(2+pw_) 
BUhelesis mt C 
an-e 2 Demon. One va 


———_, +} . (4.57) 
me (p-2){v [(p-2)]-1} 


The auxiliary parameter wv(w) is computed as follows (see Appendix-C, 
bdeive(C.18) - (C.20)): 


-1 
v(a) = [§2/w z(2)] (4.58) 


v(atqu.) for |q|>1 


=) (4.59) 
v(nw.) for |n|>1 
. . “a 
TRUE pata Ps a ae be (4.60) 
Festep SD: Thermatrjx of ‘coetficients ay -2 52 is inverted: 
fa. ] = fa 7 -M < Q S M ° (4.61) 
p-£,8 p-2 48 mae’ 


In order to obtain a matrix that relates an open-circuit noise voltage, 
rather than electric field, to the basic noise process F(w) (as in Eq. 


(3.20)), the elements a are multiplied by the corresponding relation 


P-L 52 
(see Appendix-A, Eq. (A.25)): 
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so that an expression for the normalized noise voltage eons V/V, is 


obtained: 
e j 
Weel inte) ovata de (4.63) 
where 
Vsp = v,(2tpw,) : Fe = F(2t2u) : (4.64) 


The above procedure has been safer computationally than the alternative 

of changing the original matrix before inversion. The latter alternative 
often resulted in an ill-conditioned matrix. 

Step 6: The initial spectra can now be computed from Eqs. (3.30) - (3.35), 
where, instead of Vos its normalized value ete V/VE is used, since the 


j 


matrix coefficients Chee 


are now also normalized to Vb Using Eq. (4.17) 


S p(2tew,) - el, , (4.65) 
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and choosing the dc current | (typically in the range 10 - 100 mA), we 


obtain 


: i a 
Smm(2) = 2¢me/Tg) 2 berg o/c + Sy-2,9/Me | (4.66) 
gore z€ 
S44 (2) = 2(me/1,) ey! S1-250/Me Fhe ewer Me (4.67) 
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Replacing ve IMEEGS:. «(3x3 (3.34) vand. (3.38) by Vb (ve = Vin case oF 
IMPATT diodes), we have 


Ml eas hese 
Sip (2) = 2(me/I,) jaye (4.69) 
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The spectra of the normalized open-circuit noise voltage at the sidebands, 


from Eqs. (3.36) and (3.37), are 


pe 2 
Sy(wot) = 2(me/T.) yl T-ss2! (4.72) 
Nias, 2 
Sy (4-2) = 2(me/T ) L alee : (4.73) 


The additional factor of 2 in formulas (4.66) - (4.73) has been introduced 
to obtain the magnitude of real spectra with 2 positive only. 
Calculation of the output spectra: 


In general, all the initial spectra must be inserted into Eq. (2.60). 
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However, if S54 oe ee the simplified method given in Sec. 2-5 can be 


used to compute the output spectra. The previously computed initial 


spectra Sam) » Sgt Sb (2) S..(2) and Sip (@) are then inserted into 


ob 
the equations (see Ch. II, Eq. (2.69) and (2.70) ) 
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where the coefficients pl. givenebyakass (2.63),.<(2364), (2-71) and.(2.72) ¥ 


can be written in a normalized form, realizing that 
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where psi wCyZo: 
The expression for S59 (%) can be greatly simplified by using 
Vlaardingerbroek's condition given by Eq. (2.53) which simply means that 


for IMPATT diodes 


n-v = 0 (4.81) 


1 pa 
and hence P32 = 0. 


Therefore only the primary noise remains: 


S,,(9) = [Pas |osy(a) (4.82) 


and the rms frequency deviation is 
Af = (F/O, JIS, lo) Jo (4.83) 
rms c ‘ex’ =" oo ‘ : ; 


The expression for Say (2) can be similarly simplified only if 


27 (2)>= and all modulation noise disappears. Then 
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Se) Daa ee (4.84) 


However, for Z7(2) finite, the magnitude of the modulation term in the 


be and one only has 


complete Eq. (4.74) has already been assessed elsewhere 
to be concerned about the influence of the cross-spectral term (second 
term on the right-hand side of Eq. (4.74). One might speculate that 


this term could cause a noise minimum at a bias impedance lower than 
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Z7(2)>©. However, should the cross-spectrum Sib 6%) turn out to be 
smaller than the smaller of Sam (2) and Sip (2) » such a compensation is 
impossible and only the primary noise given by Eq. (4.84) needs to be 
studied. 

Finally, for comparisons of the computed and measured spectra it 
is advantageous to plot these spectra in dependence on the diode output 
power rather than voltage since the latter quantity is difficult to 
measure with reasonable accuracy in an actual oscillator. This is carried 
out in Chapter V, Sec. 5-3. However, the voltage scale of the computed 
spectra must first be changed to output power scale, using the known 
values of series resistance Re and negative resistance Re Z(w) in Eq. 
(B.12), and using Eq. (B.9) which relates diode impedance and voltage. 
From the maximum RF power generated by the measured diode at a given dc 
current, and from the known diode impedance and series resistance at this 
power level, the corresponding rms voltage across the diode is obtained. 
This is usually close to the maximum attainable rms voltage which, 


according to several references 29#9:30 


» 1S approximately equal to 0.3V, 

for Silicon abrupt-junction diodes, while for Germanium diodes this maximum 
is approximately 0.57V, (the latter datum is from Ref. 19). The correspond- 
ing value for GaAs abrupt-junction diodes is believed to be between the 

data for Silicon and Germanium diodes. To the knowledge of this author, 

no more accurate information is available in the literature. The relative 
changes in diode impedance for corresponding changes in diode RF voltage 

can be computed from Eq. (B.9) and, assuming that the series resistance 


R. does not change with RF voltage, the power at the oscillator output can 


_ be computed for each value of the diode RF voltage. The computed and 
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measured noise spectra of the oscillator can then be graphically compared, 
using output power as the x-coordinate and noise spectral density as the 


y-coordinate. 


4-3-2 Plots of Computed Spectra and Discussion 


The procedure described in the previous section yielded data for 
the initial and output noise spectra as a function of the physical and 
electrical parameters of the IMPATT diode and of the external oscillator 
circuit. Numerous plots can be constructed using the above mentioned 
data. The most important of these are shown in Figs. 4-2 to 4-9 in this 
section. With the exception of Figs. 4-8 and 4-9, which compare the 
output amplitude and frequency noise spectra of IMPATT diodes made from 
Si, Ge and GaAs, these graphs are for Si diodes and depict the initial 
auto- and cross spectra, or the open-circuit noise voltages of these, 
in dependence on diode and circuit parameters. These are: RF electric 
field, RF voltage, bias circuit conductance, drift transit angle, relative 
avalanche region width and avalanche-to-carrier frequency ratio. 
Corresponding plots for GaAs and Ge diodes are similar and need not be 
shown here. 

The frequency dependence of the computed initial spectra turned out 
to be negligible in the baseband frequency range of interest. Hence, the 
parateter Of 0, is not included in most of the plots. All noise spectra 
have been computed for a de bias current of 10 mA. For other currents, 
the power spectra and the rms frequency deviation must be divided by 


Wz 


I (mA)/10 and by [1 (mA) /10] » respectively. 


The range of the diode parameter w./Wq (avalanche-to-drift region 
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length ratio) has been chosen somewhat arbitrarily around Kuvas' values 
in Table 4-1. There is no need for great accuracy in this respect, 
since these values vary from diode to diode, being dependent on diode 
structure and background doping concentration“®, Moreover, the above 
ratio seems to depend strongly on how the equivalent avalanche region width 
is defined. For example, the data in Ref. 27 differ completely from those 
in Ref. 28, where a different definition of avalanche region width is used. 
Yet, these uncertainties do not detract from the value of the plots given 
here, since the plots also show how the noise spectra are affected by the 
changes of the parameter involved. 

The drift transit angle Wot is also affected by the uncertainty 
in avalanche region width; it further depends on the oscillation frequency 
and may vary from below 2 to more than 3 for some structures. The value 
of 3 has been chosen for some of the plots with a fixed WoT since, 
~ according to Ref. 20, this value should be an optimum compromise between 
good noise performance and reasonable output power. 

Fig. 4-2a shows the dependence of the open-circuit noise voltage 
spectrum on the parameter (m/w t)E.,/E, for several values of w/w. 
While decreasing w/a results in decreasing noise voltage spectral density, 
it also makes the dependence of the latter on electric field steeper. 
This can be better seen in Fig. 4-2b, where the open circuit noise voltage 
spectrum, normalized to its small-signal value, is again plotted as a 
function of the quantity (m/w) 2 Ea eD with wh fue as a parameter, 

20 


and a comparison is made for we/ue = 0.6 with the curve of Kuvas -. The 


increasing difference of the two corresponding curves at higher fields may 
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be caused by the different matrix sizes used (M=5, compared to M=2 of 
Kuvas). It was found in our calculations that a value of M larger than 
5 did not influence the results significantly. Also, there is some 


a that Kuvads' noise voltage increases too rapidly 


experimental evidence 
with increasing electric field. However, its basic behaviour is the same 

as that of the present data in Fig. 4-2. At low RF fields, the open- 
circuit noise voltage does not differ appreciably from that existing under 
pure dc excitation; at higher fields, it increases rapidly with increasing 
field. 

Fig. 4-3 shows the dependence of the initial modulation spectra on 
the normalized RF diode voltage, for various values of wh /ue aS a parameter. 
At very low RF voltage, the amplitude-and phase spectra are almost equal; 
to this, there corresponds an open-circuit noise voltage in Fig. 4-2 
almost equal to that existing under pure dc excitation. With increasing 
RF voltage (and electric field) the difference betwen the initial amp 1itude- 
and phase spectra increases also. Figs. 4-2 and 4-3 thus confirm the 
general considerations in Chapter II; the perturbing noise voltage source 
is amplitude-dependent and its orthogonal components v(t) and ve(t)s which 
are directly related to the initial modulation rates have, in general, non- 
equal spectra. Also shown in Fig. 4-3 is the RF voltage dependence of the 


cross-spectrum of amplitude and phase, S_.(2), in absolute value. It is 


md 
seen that, for a high impedance in the bias circuit, Z(2)>~,this cross- 
spectrum is insignificant. 

In Fig. 4-4, the normalized baseband spectrum Sip (2) and the absolute 


value of the amplitude-baseband cross spectrum Sab (2) are plotted as 
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functions of the normalized RF diode voltage Vi/Vp 3 for infinite impedance 
in the bias circuit. Two inter-related phenomena are indicated: 1) even 
for infinite bias circuit impedance there exists some correlation between 
amplitude fluctuations and baseband noise voltage, and this correlation 
increases with increasing RF diode voltage, and 2) the baseband noise 
voltage spectral density also increases with the latter. 

Figs. 4-5 and 4-6 show the dependence of all the initial noise 
spectra of interest on bias circuit conductance, at a high RF diode 
voltage. The conductances to the right of the vertical dashed line are 
not physically realizable because of the thermal resistance of IMPATT 
diodes. The position of this line has been derived from the data on 
thermal resistance and chipcapacitance of typical Si IMPATT diodes 
designed for operation at 6 - 12 GHz. It is an average position and there 
may be deviations in both directions for each individual diode. The minima 
of initial AM noise, Sam 6) 2 in Fig. 4-5 cannot be observed in practice 
since, at the corresponding values of bias circuit conductance, the 
"modulation noise", given by the third term on the right-hand side of Eq. 
(4.74), completely dominates the oscillator output. This is easily seen 
by inspecting the magnitude of the coefficients Pr, and ae io Eqm. 4.74) 
and that of the noise spectra Sim? Spb and Smab in Figs. 4-5 and 4-6. Hence 
it is concluded that the output amplitude noise is minimum at Z(2) = ~, as 
expected. 

As mentioned in Sec. 4-3-1, the output phase-or frequency noise is 
not influenced by modulation. The slight increase in FM noise with 


: : : : A : ; 16 
increasing bias circuit conductance, which is sometimes observed ~, may 
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be caused by the increase in the initial spectrum S,,(2) indicated in 


o> 
Fig. 4-5. 

Figs. 4-7 shows the influence of the diode drift transit angle and 
that of the avalanche-to-drift region ratio on the initial amplitude- and 
phase noise spectra, for three values of the RF diode voltage. An 
optimum range of WoT qs 4.5 to 5.5 radians, is indicated in Fig. 4-7a for 
Sam (2) and S454 
rapid increase in noise, especially at large RF voltages. Note also that 


(2). A deviation from this optimum range results in a 


WoT influences the relative magnitudes of the two initial spectra. The 


above optimum range of 4.5 - 5.5 radians is not in agreement with the opti- 
mum value of 3.6 radians resulting from the power/noise considerations 


21 who, from a 


of Kuvas@?, but agrees well with the optimum of Goedbloed 
small-signal theory,derived a curve similar to those in Fig. 4-7a. 
However, an actual diode operated at a transit angle in the range 4.5 - 
-5.5 radians will generate low RF power, as the optimum of the latter is 
approximately at Wot] = 2.3. From this additional consideration, the 
value of uty=3.6, given by Kuvas“?, will thus likely be an acceptable 
compromise. | 
The dependence of the initial amplitude- and phase noise spectra 
on the ratio w,/Wq is more straightforward. It is clearly seen that a long 
avalanche region is desirable for low noise, as both amplitude and phase 
noise spectra increase rapidly at small values of w/Wq- However, the 
validity of this conclusion is limited somewhat because of the Read model used. 
Finally, in Figs. 4-8 and 4-9 a comparison is given between the noise 
performance of Si, Ge and GaAs IMPATT diodes. For this comparison, the output 


amplitude noise spectrum Sum (2) and the rms frequency deviation ae have 
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been chosen. Unlike the power/noise comparison of Kuvas“?, lite TSu NOt 
attempted here to relate noise performance to attainable power. The 
emphasis is on an absolute comparison of noise modulation rates. Also, 

for a fixed transit angle, the effect of the relative avalanche region 
width is shown. For Ge and GaAs, the voltage scale differs from that used 
for Si, since saturation in the former occurs at higher RF voltages. It 

is seen clearly that information on the noise performance of diodes made 
from different materials, based on small-signal noise data or theory, may 
be rather unreliable when applied to a high-level oscillator. The effect 
of the RF voltage and that of the ratio of w/Wq is relatively small for 

Ge diodes compared to Si and GaAs. Other than that, no generally valid 
comparison can be made at high RF voltages unless it is understood that the 
diodes will have certain ‘typical’ ratios of W/Wq: For the ratios of 
w,/Wq given in Ref. 20, Figs.4-8 and4-9 show that at low RF voltages 
Silicon diodes exhibit a higher noise than Ge and GaAs diodes, in that 
order; at very high RF voltages, a GaAs diode may exhibit higher noise than 
a Si diode. A typical GaAs diode is seen to exhibit the lowest noise of 
all three, if it can be operated at a relatively low RF voltage near the 
noise minimum. 

In order to verify the theory presentedin this thesis, the computed 
output noise spectra must be compared with some measured data of actual 
IMPATT diode oscillators. For this purpose, no "typical" parameter values 
can be used. The diodes tested must be examined for these parameters in 
each specific case so that proper values are used in the computations. 

The computed spectra obtained in this way are presented, together with the 


measured data, in Chapter V, Sec. 5-3. 
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CHAPTER V 
NOISE MEASUREMENTS ON IMPATT OSCILLATORS 


The emphasis in this project has been on developing a large-signal 
theory of oscillator noise for a certain class of negative-resistance 
oscillators. This theory has been applied to IMPATT diode oscillators 
in order to demonstrate how the general theory can be adapted to a 
particular device and, at the same time, to obtain specific data which 
might explain some of the large-signal noise effects observed in IMPATTS 
or, which may even suggest the existence of some noise phenomena not yet 
observed. 

In this chapter, specific computational results will be compared 
with the pertinent experimental data on IMPATT oscillator noise. The 


data of Goedbloed*| 


have been supplemented by additional measurements 

_ carried out by the present author. For this purpose, a sensitive noise- 
measuring system had to be assembled and thoroughly tested fowiits 
limiting sensitivity (threshold) and accuracy. A direct-detection noise- 
measuring system has been chosen. It is described briefly in Section 5-1. 
The most important data of the IMPATT oscillator and of the diodes used 
are summarized in Section 5-2, while Section 5-3 contains both the 


pertinent theoretical results of this study and the available experimental 


evidence. 


5-1 The Noise-Measuring System 


5-1-1 Small-Signal Measuring Principles 


Direct-detection noise-measuring microwave systems have become wel] 
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established lately because of such advantages as better sensitivity (if 
sufficient input power is available) and the avoidance of complications 
arising from the use of a local oscillator in a superheterodyne system. 
A spectrum analysis is obtained separately for AM and for FM noise 
modulation, so that the AM detector is required to not respond to FM, and 
vice versa. 

A good discussion of cman iaetgnal measuring principles is found in 
Refs.32 and 33. Basically, the waveform is assumed to be expressible in 


the form 
v(t) # veel + m(t) Jcos[w, + Aw(t)]t (5.1) 


for Aw(t) small and changing slowly with time (quasi-stationarity). An AM 
detector will yield an ac output dependent only on m(t) and will not 

_ respond to changes in frequency Aw if its bandwidth is sufficiently large 
compared to peaks of Aw. Similarly, an FM detector normally Wreids an 


output voltage 
v(t) = K-V [1 + m(t)]"aw(t) , n= lor 2 (5.2) 


and thus v will be a measure of Aw if m(t) << 1. This condition is usually 
well satisfied when m(t) is noise modulation only. 

The basic direct-detection measuring system described here, 
analogous to those of Refs.32 and 33, is shown in Fig. 5-1. Its operation 


is as follows. 
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When channel II is interrupted by the switch SW-1 and the E arm 
of the hybrid tee match-terminated, the hybrid tee and diodes D, and D, 
function as an ordinary AM detector (either square-law or linear), if the 
diodes are connected in phase. If they are connected in opposition, the 
two diode signals will cancel if the two detector arms are identical 
electrically. The degree to which this cancellation is achieved is 
expressed by the suppression factor (SF) 


SF(dB) = 20 log( (ons) 


y inbaw! (Ban ; 


where V is the output voltage in unbalanced condition, and Hear is 


unbal 
the output voltage in balanced condition (diodes in opposition). Diode 
noise that is uncorrelated between D, and D., will not cancel and wili 
determine, together with the amplifier noise, the AM threshold of the 
system, if the suppression factor is sufficiently high. 

For FM measurements, both channels of the system shown in Fig. 5-1 
are used. The carrier suppression filter consists of a circulator and a 
high-Q reaction cavity. Ideally, the carrier is completely suppressed, 
which requires that when the cavity is tuned to Wes it is exactly 
critically coupled. If the frequency deviates from We by a small amount 
Aw, a Signal proportional to Q,(Aw/w.)s shifted in relative phase by 
£2908, will be transmitted to the E arm of the hybrid tee; this signal will 
be much smaller than the reference signal (containing both carrier and 


sidebands), transmitted to the H arm of the hybrid tee. With the diodes 


switched in opposition, this circuit operates as a phase-sensitive 


ae 


7 


mis 3 sag bre [We dadkwa etd xd oa at 


at pg 4M bis /g tt rbnoo beons fsdew th sot ugha sti af fedm” 

shod .(notttzeqgo Wt esbots) nord this’ badnstsd 9 nt opatfov uate 

[Tiwebes Toons. ton rte <A bis (0 nsowded ‘bate feryeony af reas 

ont to blomzenrie MA. silt -a2ton bibs a ame : 
A eet Pane totttwe 2t sodas} note 1qqu ian tie oil 

[-2).017 it mide mane (anit: pene tot a momen wasem MA 04 me 

6 bas Yors hurt & a edettanoo mot? notee re . vs * 69 oni 2 ° 

a bseesyqque ad a erly a, 5 hy 2 

yisosxe 2k If of. bent at Wi Fea, tt pe ted? bh ns * 

Jmwoms Themes Na yo MoH aate tab Yonsupert ent #1. .betquos ylfsakstao | 
ve s24dQ° 9V tsetonint hegre. «(ows gp: % tenot sroqong, fsngte oe 
(tw. Tenpte zig poe ey ‘og tome 2 old CT ” ; | oo aay 


fourm A od 4 


ham 


ashore a AsiH 298) gt: ad 40. ints Wa ‘be 


~ 


oviti renge~ -926da, 5 28 283 s%9q0 stuoys 2a 


oy 


detector, producing an output voltage proportional to (amplitude in 

channel II) cos(>, + 1/2) (for linear detector), or to (amplitude in 
channel I) (amplitude in channel IT) cos($., + 1/2) (for square-law 
detector), where by is the phase difference between the two channels. If 
o, is made equal to 1/2 by means of a phase shifter in channel I, we obtain 


maximum output voltage 
v(t) = K[1 + m(t)]" + Qoaw/o (5.4) 


where n = 1,2 for linear and square-law detectors, respectively. Thus 
Vib) iaohuhei fan(t)o<<.1. 

Once the system has been calibrated, its FM threshold, i.e., the 
lowest measurable Aw, can be determined from the threshold voltage obtained 
previously with channel II interrupted. 

The detector output voltage, proportional to m(t) or Aw(t), is 
“amplified and subjected to spectrum analysis. The root-mean-square value 
of m(t) or Aw(t) in a given frequency slot (usually 100 Hz or 1 kHz) in 
the range of modulation frequencies of interest is then measured. For 
maximum system sensitivity, the spurious noise must be kept at a minimum, 
since the threshold voltage is determined by the noise added by the 
detector diodes, the subsequent amplifier, and possible spurious voltages 


due to ground loops and other extraneous causes. 


5-1-2 Threshold and Accuracy 
When checking the AM threshold of the system in Fig. 5-1, an 


unexpected rise in threshold below approximately 10 kHz may be encountered. 
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As mentioned in Ref. 32, such a rise in threshold in superheterodyne 
systems is caused in some cases by excessive FM of either the signal 
under test or of the local oscillator, in conjunction with FM-to-AM 
conversion in the IF amplifier. A similar phenomenon can occur in the 
case of a direct-detection system if the FM noise of the oscillator under 
test is relatively large and the detector used is narrow band. This is 
the case with some envelope detectors because of their comparatively high 
RC, product. 

A series of experiments has been carried out by this author, which 
have shown that the threshold at low frequencies can be reduced by broad- 
banding the detectors. This procedure, together with a more detailed 
description of the measuring system, is explained in Ref. 34. Since these 
problems are auxiliary to the main theme of this thesis, only the results 
will be cited here: with a detector Q of about 20, and an input RF power 

of some 40 mW, an AM threshold range of -150 to -155 dB/100 Hz has been 
; attained at modulation frequencies from 100 Hz to 50 kHz. With an input 
power of 15 mW and with the detector frequency response adjusted flat 
within 1 dB over a 100 MHz range, the AM threshold is approximately 
-140 dB/100 Hz (at 100 Hz) to -155 dB/100 Hz (at 50 kHz). This is 
sufficient for noise measurements on IMPATT diode oscillators. 

Another problem associated with systems of this type is the effect 
of their various imperfections on measurement accuracy. The operation of 
this ciruit and its calibration have so far been described in terms of 
assumed "ideal" conditions; namely, the hybrid tee is perfect, the cavity 
is tuned precisely to the carrier frequency of the RF signal, and the 


cavity input impedance is perfectly matched at that frequency, thus 
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producing no reflection of the carrier. In practical measurement setups, 
conditions differ significantly from the ideal situation causing various 
AM-to-FM and FM-to-AM conversions. A computer study of these effects has 
been carried out by Fikart et Hes for certain typical values of 
oscillator AM and FM noise and for a certain typical range of deviations 
from the above mentioned ideal conditions. This study is explained in 
more detail in Ref. 35. The results, shown in Fig. 5-2, indicated that 
FM noise measurements can be made accurately for modulation frequencies 
within the reflection cavity bandwidth, whereas AM measurements using the 
"“double-channel mode" are highly unreliable. Hence for AM measurements 


the single-channel mode must be used only. 


5-1-3 Complete Measuring System 


The complete system for noise measurements at 6 GHz is shown in 
Fig. 5-3. It is essentially the same system as that in Ref. 33 except 
that it uses Schottky-barrier diodes as envelope detectors, as in Ref. 32, 
instead of small-signal backward-diode detectors. Some instruments have 
been added for convenience: a 4-8 GHz sweeper, an oscilloscope to permit 
fast adjustment and retuning of the detectors, and a frequency counter for 
direct frequency observation. The video detector output has been arranged 
as shown in Fig. 5-4. The two diodes are reversed with respect to each 
other in their mounts and when threshold measurements are taken (or FM 
operation required) they are simply connected in a bridge configuration. 
The two variable resistors facilitate balancing of the bridge since, 


generally, the detector video impedances are not equal when the voltages 
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are, and vice versa. Their value is chosen so as not to greatly 
influence the total noise. 

For actual AM measurements only one detector is calibrated and 
used. This procedure has the advantage of removing spurious low- 
frequency voltages due to an ungrounded detector mount (when the detectors 
operate in parallel or in series). Hence, total battery operation is not 
then required; only the low-noise amplifier is battery operated. 

For this system, operating at 6 GHz, a very low-noise Schottky- 
barrier diode pair HP 5082-2354 has been chosen. The amplifier used has 
an equivalent rms input-noise voltage of about 0.02 - 0.03 uV in a 100 Hz 
band at low frequencies. 

During the measurements, the system has been operated at an input 
power of 15 mW and with the detectors tuned for a flat response in a 
100 MHz range around 6 GHz. The AM threshold was from -140 dB/100 Hz (at 

modulation frequency 100 Hz) to -155 dB/100 Hz (at 50 kHz). The 
corresponding FM threshold, using a reflection cavity with Q = 15000, was 


approximately 0.02 - 0.03 Hz/100 Hz (at 50 kHz). 


5-2 Auxiliary Data: IMPATT Diodes, Coaxial Resonator and General 


Measurement Procedure Used. 


Commercial Si and GaAs IMPATT diodes were used in the noise 
measurements. After completion of the noise measurements, the physical 
dimensions of the diodes were examined by measuring the semiconductor 
chip area, and the diode capacitance dependence on reverse voltage. The 


determination of chip area turned out to be impossible in the case of 
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the GaAs diodes because of their special contacts, so the data from 
Ref. 28 were used to calculate the depletion layer width from the 
known breakdown voltage. The combined physical and electrical parameters 


of the two types of IMPATT diodes at 6 GHz are given in Table 5-1. 


TABLE 5-1 PARAMETERS OF MEASURED IMPATT DIODES AT 6 GHz 


The diodes were inserted into a tunable coaxial resonator, a sketch 
_ of which is shown in Fig. 5-5. The two movable slugs made possible fine 
frequency tuning and variable load coupling. The diode holder is of 
standard design (see, e.g., Ref. 37). Its primary function is that of a 
good electrical contact with the diode and adequate heat removal. 

The diodes were biased bya high-impedance (current) source, typically 
in the range of bias currents 20 to 30 mA. At a given dc current, the RF 
power was adjusted by means of the tuning slugs to several different 
values. The output signal was monitored on a spectrum analyzer to ensure 
that no spurious oscillations were taking place (these cause a large 
increase in both AM-and FM noise, compared to single-frequency oscillation). 

For each output power adjustment, the external Q of the resonator 


was measured by the injection locking method?®, using an auxiliary 
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generator and circulator. After that, the IMPATT oscillator was connected 
to the noise measuring system and AM-and FM noise measurements taken. The 
baseband noise voltage across the bias leads was also measured. During 
these measurements, the oscillation frequency was kept constant, ce = 5.94 
GHz. The detected AM-and FM noise spectra, as well as the baseband noise 
voltage spectrum, were measured in the rangeof baseband frequencies 
100 Hz - 50 kHz. This frequency range enables one to determine if any 
flicker noise is superimposed on the white noise spectrum, normally 
expected from IMPATT diodes. The measured spectra are compared with the 
pertinent computed spectra in the next section. 

In the course of measurements, the two GaAs diodes used were found 
to be somewhat unreliable compared to the Silicon ones. One was damaged 
at an early stage, apparently by a pulse in the bias circuit, when baseband 
noise voltage was being measured. Although this particular measurement 
was made repeatedly in the case of the Si diodes, the latter did not suffer 
any noticable damage or deterioration. The other GaAs diode, on which 
baseband noise measurements were not made, lasted through the measurements, 
but its performance gradually deteriorated, without any apparent external 
cause. This aging manifested itself mainly in a decreased maximum 
available RF power at the end of the measurements. On examining this 
diode, a change in its breakdown voltage was found, from the original value 
of 80 volts to 60 volts. Also, a relatively poor contact between the 
semiconductor chip and the metallic stud was discovered; this condition 
normally causes an excessive junction temperature and contributes to the 


generation of excess thermal noise! This condition is believed to be 
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partly responsible for the unusually large FM noise measured in this case. 

Normally, GaAs diodes, if well manufactured, are known to exhibit quite 

low FM noise?” at not too low baseband frequencies. In any case, the 

above-mentioned deficiencies appear to be more often encountered in 

GaAs diodes than in Si diodes, the technology of which is highly developed. 
In these circumstances, a comparison of the theoretical- and measured 

dependence of noise on output power is not representative. For this reason, 

only the baseband frequency dependence of both AM-and FM noise of the 

other GaAs diode will be shown, in comparison with the corresponding data 


of the Si diode measured. 


9-3 Comparison of the Computed Noise Spectra with Experimental Data 


One of the interesting computational results in Section 4-3-2 is the 
_ generally different rate of change of the amplitude-and frequency noise 
spectra with a change in diode RF voltage. As mentioned in Chapter ee if 
Kuvas' large-signal analysis were used@9, the two noise spectra would 
exhibit an equal dependence on diode RF voltage. The experimental results 


of Goedbloed | 


show clearly that, while FM noise increases relatively 
rapidly in the high voltage region, the AM noise there almost does not 
increase at all. Fig. 5-6 shows Goedbloed's results in comparison with 
the theoretical curves obtained by the present author. Also shown in 
Fig. 5-6a is a theoretical curve by Goedbloed (small-signal theory). 


The comparison is not absolute since not all the diode and circuit 


parameters can be reliably derived from the data given in Ref. 21. 
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Fig. 5-6a FM NOISE AS A FUNCTION OF RF VOLTAGE 
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Note that the data in Fig. 5-6 were obtained with a constant load 
resistance, i.e., the RF voltage across the diode was varied by varying 
the bias current. 

In Fig. 5-7, the influence of drift transit angle on FM noise 
is shown, as predicted and measured by Goedbloed@! , and as computed from 
the large-signal theory of the present author. As before, the comparison 
is not an absolute one, because of some missing data in Ref. 21. However, 
the position of the minimum at 4.5 < WoT gs 5.5 seems to be little dependent 
on other diode and circuit parameters. This conclusion can be drawn from 
the results of Sec. 4-3-2 of this thesis, from Goedbloed's additional 
observations¢| and, also, from the data in Ref. 40. 

Figs. 5-8 to 5-11 summarize the experimental results of this author, 
in comparison with the pertinent data obtained from the computations 
described in Chapter IV. In these figures, the noise spectra are related 
to output power, a more useful parameter, rather than to diode RF voltage. 
| Fig. 5-8 shows the theoretical and measured dependence of AM-and FM 
noise on output power. An important difference of these experimental 
results with respect to Fig. 5-6 is the tuning procedure used. The data 
in Fig. 5-8 were taken at a constant bias current, but with variable load 
resistance. This was achieved by adjusting the tuning slugs to obtain 
the required RF power, while keeping the frequency constant. It was found 
that values of relative output power smaller than approximately 0.65 were 
difficult to obtain without the presence of spurious oscillation. Further, 
an effect similar to that seen in Fig. 5-6a was observed; namely FM noise 


was larger than expected from theory at low power levels, with spurious 
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Fig. 5-7 FM NOISE AS A FUNCTION OF DRIFT TRANSIT ANGLE 
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Fig. 5-8a FM NOISE AS A FUNCTION OF OUTPUT POWER 
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oscillations present. It is possible that the experimental data in 

Fig. 5-6a were obtained under such condition, since spurious oscillations 
may, without special attention to them,easily go unnoticed. Fig. 5-8a 
shows that the measured FM noise increase in the high-output power region 
is well described by the theoretical curve. The theoretical curve of 

AM noise exhibits a little steeper rate of increase in the high-power 
region, than the measured deen On the whole, the agreement between 
theory and measurement is very good. 

The measurements described above were also taken for bias currents 
of 20 and 30 mA; the results obtained were essentially the same as those 
in Fig. 5-8. This is illustrated in Fig. 5-9, which shows the FM-and AM 
noise dependence on bias current, both theoretical and experimental. It 
must be stressed once again that for each bias current the oscillator was 
tuned to the same relative output power. The resulting dependence of AM- 
and FM noise on bias current is very flat, unlike that observed in the 
ae: of a constant load resistance, with variable output power. In the 
latter case, the FM noise has a distinct minimum for a certain bias current 
(this can be derived, e.g., from Fig. 5-6a). 

The baseband noise voltage spectrum is normally of little interest 
in oscillator noise theory; nevertheless it was measured in this case in 
order to determine whether the down-conversion predicted from the theory 
in Sec. 3-4 and indicated in Fig. 4-4, actually takes place. As is seen 
in Fig. 5-10, this is indeed so, although the theoretical noise voltage 
spectrum, compared to the measured one, increases more rapidly with RF 


Output power. 
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Finally, Fig. 5-11 shows how "white" the noise spectra of IMPATT 
diodes are. Computations based on the present theory resulted in flat 
noise spectra in the range of baseband frequencies shown in Fig. 5-11. 
As expected, the Si diodes measured exhibited noise spectra which were 
effectively flat to a frequency as low as 100 Hz. The GaAs diode 
produced AM-and FM noise spectra which increased gradually at frequencies 
lower than approximately 10-20 kHz. This agrees with Ref. 41. This 
"flicker noise" phenomenon is known to be associated with the quality of 
the semiconductor surface. At present, a surface quality comparable to 
that of Si diodes is difficult to achieve in the case of GaAs devices. 
Therefore, any theory of noise for the latter devices that does not take 
into account flicker noise, can hold only for baseband frequencies 


greater than approximately 20 kHz. 
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CHAPTER VI 
SUMMARY AND CONCLUSIONS 


In this thesis, a large-signal theory of noise in oscillators 
has been presented, using a modified quasistationary perturbation theory 
and an RF-amplitude dependent perturbing noise voltage (current) source. 
The fluctuations in amplitude ste phase are assumed RF-amplitude 
dependent, but still sufficiently small to permit making the usual 
quasilinear approximation in the perturbation theory of oscillator noise. 
The method introduces the concept of initial modulation rates of AM-and 
FM noise, whose auto- and cross spectra can be determined from a differ- 
ential equation which described the active element via open-circuit noise 
voltage or short-circuit noise current at the sidebands about the 
oscillation frequency. The spectrum of the baseband noise voltage (current) 
can be also determined in this way, as well as the cross spectra of the 
latter with the noise voltage (current) sidebands. The above mentioned 
auto- and cross spectra depend on the magnitude and phase of the periodic 
voltages and currents present in the circuit as well as on the circuit 
impedance at baseband frequencies. The AM-and FM noise spectra at the 
oscillator output are computed from the initial spectra by means of a 
modified quasistationary perturbation theory. Unlike the large-signal 


20. the present approach makes it possible to separate the 


theory of Kuvas 
computation of the AM-and the FM noise. This is important, since it turns 
out that the dependence of AM noise on signal level in oscillators is 


generally different from that of FM noise. The method is suitable for 
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active elements that can be described by a differential equation in 
conjunction with a Langevin noise source. 
The theory has been applied to IMPATT diode oscillators, using the 


20 This reference has 


Read diode equations in the form given by Kuvas 
also been used for the calculation of several constants needed to 
approximate the behaviour of abrupt-junction diodes made from Si, Ge and 
GaAs. Both the initial and output AM- and FM noise spectra have been 
computed, for various combinations of IMPATT diode and circuit parameters. 

A noise-measuring system at 6 GHz has been assembled and tested. 
Two types of IMPATT diodes were operated as oscillators in a tunable 
coaxial resonator. The AM- and FM noise spectra of these oscillators 
were measured by means of the above mentioned system, for different bias 
currents and output power adjustments. 

The computed spectra were compared with the experimental data of 


Goedbloed<! 


and those obtained by the present author. A good agreement 
between the theoretical and experimental data has been found. “In addition, 
the following noise properties of IMPATT diode oscillators have been 
revealed or confirmed: 

1) The perturbing noise voltage is strongly RF amplitude dependent; its 
two components, one in-phase and the other in quadrature with the oscillator 
signal, may, in general, be correlated and have non-equal spectra. The 
concept of an independent perturbing noise voltage source is valid only at 
small RF voltage amplitudes. 

2) Minimum AM-and FM noise at the oscillator output is obtained for a 
large impedance in the bias circuit. 

3) Both FM- and AM output noise modulation spectra first decrease and then 
increase again with increasing RF diode voltage at the frequency of 


oscillation. Typically, the minimum of the AM noise occurs at a higher 
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RF voltage amplitude than that of the FM noise, and the subsequent increase 
in the AM noise is less rapid than that of the FM noise. The rate of 
increase, positions of the minima and the absolute values of AM- and FM 
noise spectra depend on the diode material and on the physical and 
electrical parameters of the diode. Experiments have shown that, for 
typical commercial Si IMPATT diodes, a decrease of some 30% from the 
maximum RF power at a given bias current results in a considerable 
reduction of the FM noise, approximately by 10 dB, while the AM noise 
typically changes by only a few dB. 

4) A low ratio of W/o: large drift transit angle up to approximately 

5 radians and a wide avalanche region relative to the drift region reduce 
both AM-and FM noise. The latter increase very rapidly with decreasing 
avalanche region width at high RF voltages. For minimum AM-and FM noise, 
a given diode should be operated at a high oscillation frequency. This 
decreases the ratio w/o, and increases drift transit angle. 

5) Comparison of the noise performance of oscillators using IMPATT diodes 
made from different materials (typically Si, Ge and GaAs) is difficult to 
carry out in general terms. Small-signal noise theories indicate lowest 
noise in GaAs diodes, with Ge and Si diodes exhibiting higher noise, in 
that order. The theory presented in this thesis shows that these 
conclusions are not necessarily valid under large-signal conditions, where 
small differences in diode dimensions may have a greater influence than 
the material itself. For typical IMPATT diode dimensions such as those 
given in Ref. 20, the theory indicates that, at low RF voltages, GaAs 
diodes have a lower noise than Ge and Si diodes, in that order; at very 


high RF voltages, a GaAs diode may exhibit even higher noise than a Si 
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diode. Finally, the potentially good properties of GaAs diodes at medium | 
RF voltages are somewhat compromised by their relatively high flicker 
noise at baseband frequencies below approximately 20 kHz. 

It is felt that the above results and the corresponding plots are 
useful as guidelines for the choice of IMPATT diodes and their operating 
conditions. They may be useful especially in those applications where 
suppression of FM noise by simple "internal" means would be sufficient, so 
that such ‘external " means of noise suppression as injection locking and 
high Q cavity stabilization need not be used. In this respect, the 
particular results of this study complement the power/noise considerations 


of Kuvas<o 


» while the general theory presented in this thesis makes it 
possible to estimate large-signal AM-and FM noise of some other oscillators 


as well. 
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DERIVATION OF SOME READ DIODE RELATIONS IN NORMALIZED FORM 


The purpose of this Appendix is to derive the relations between 


electric field and external current, conduction current and terminal 


voltage for the Read diode, in a suitably normalized form. 


We shall start with the introduction of the following normalized 


quantities in the frequency domain: 


Normalized conduction current 


i (w) = I (w)/1, 


Normalized external current 


ig, (w) = I, (w)/1, 


Normalized electric field 


Normalized voltage 


V(w) = V(w)/V 


(A.4) 


In the above, the frequency dependent quantities, denoted by capital 


letters, are either Fourier transforms or Fourier series components of 


the corresponding stochastic or deterministic quantities, respectively, 


in the time domain. 


Furthermore, we shall use the "avalanche frequency" parameter 
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W. defined by 


eae 
a 


Ww mI ,/(eAtE,) = mI ,/(eAkt .E, ) 7 (A.5) 
where I, is the dc (bias) current through the diode, E, is the electric 
field at breakdown, « is dielectric permittivity, A is the effective 
cross-sectional area of the diode, e. is the transit time through the 
avalanche region, k is the correction factor as explained in Sec. 4-1-3 
and m is the exponent in the approximation of a(E.) in Sec. 4-1. 


A critical parameter for Read diode oscillators is the ratio of 


the avalanche-to-carrier frequency: 


(w/o) = ml / (eAke E06) = (m/2k) (1/u.t,) (15/0 CV) 
(D/2)(Tj/u CyVy) (A.6) 


~ where Cy is the diode geometrical capacitance 
Cy = cA/W , (A.7) 


Vi is the breakdown voltage of the abrupt p-n junction which jis to be 


approximated by the Read model (c.f. Sec. 4-1-3, Eq. (4.25)) 


Vi = Epw/2 (A.8) 


and the parameter D is given by 


D=m/(kwt,) . (A.9) 
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D can be more conveniently expressed in terms of other important Read 
diode parameters; the drift transit angle Wotg where Tq is the transit 
time of electrons or holes through the drift region, and the equivalent 


avalanche-to-drift region length ratio wa/Wa: 
D = (m/k)/(woty x w/w) : (A.10) 


where Wa and Wy are the avalanche- and drift region length, respectively. 


Let us now turn to Eq. (4.12) for the external current: 


dE. Ge) 
(tine de(t \ytyeA.— 


é Taian (A.11) 


This relation is valid at the boundary of the avalanche and drift regions 


and is a special case of the more general equation 
: dE (t,x) : 
(t) = I (t,x) + cA dt : (A.12) 


Eq. (A.12) can be integrated along the drift region: 


W W W 
1 Z elt eA d 
a | Loy (t)dK = Toy(t) = Gf Teltaxdde + Gael E(tax)dx) 
a Wa na (A.13) 


Now the voltage V(t) across the diode is given by 


W W 
v(t) = - feCt onder = -£,(t)w, my Fipscdx & (A.14) 
(@) W 
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since the electric field may be assumed to be independent of position 
in the avalanche region. 

Because of the absence of ionization in the drift region and with 
the assumption of a saturated drift velocity of the carriers, the particle 


current I (tx) can be expressed as a travelling wave 

i At.) )= I t-(x-w.)/v, ox = wd ! (A.15) 
With a new variable 
: (A.16) 


the first integral on the right-hand side of Eq. (A.3) can be rewritten 


as follows: 
: W V. tery UE 
— |] I_(t,x)dx = - — I (t')dt' = — Patol ate RZ) 
Ne iE t-ty 
where I(t’) is the conduction current at x = Wo 
Eq. (A.13) now becomes 
1E 
Gye tet 1 (tide ek ee Vt) e E(t) (A.18) 
ex Tq Cc d Wy dt a a : F 
t-ty 


Let us Fourier transform Eqs. (A.11) and (A.18): 


(w) = I.(w) + jweAE,(w)  . (A.19) 
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Normalizing and using Eq. (A.6) we obtain 


i__(w) = i (w) + J (weAE,/T, )e, (o) = 


_ 7 + gc 
= i (w) + §(w/o,)D(w,/w,)“es(w) (A.20) 


From Eq. (A.18) we have 


1-exp(-jurt 4) 


ON faces bet Sot g a AVC) HES (aie 


= @(w)1(a)-d(w/v Ju Cylw/wy)V(w) + E(w)w,J 9, (A.21) 


where 


A(w) = (1-exp(-juty)/juty : (Ae22) 
Normalizing and using Eq. (A.6) again, we obtain 


jax (0) = 8(a)i (w)-5 (w/wg) (w Cy/Ty) (W/Wg) (E,W/2)v(w) + 


- J (w/w Jw Cy/T,) (WE, ) (wa/wyle, (w) 


= @(w)Iq(w) - §(w/w,)(D/2)(w,/w,) “(w/wg)v/u) + 


C 


~ j(w/w,)D(w,/o,)“(wa/wyleg(w) (A.23) 


where 6(w) can also be written in the normalized form 
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6(w) = (1-exp[-j (w/w Jo tg] Li (w/a Jorg] ; (A.24) 


On eliminating i. (w) from Eqs. (A.20) and (A.23) and setting 
ig, (w) = 0, an equation is obtained for the normalized open-circuit 
voltage v (w) in terms of the normalized field e,(w): 


V"(w) = -2Le(w) + (wa/mg)ILI+(Wa/wg)]le,(o) (A.25) 


Similarly, eliminating ig, (o) from Eqs. (A.20) and (A.23) yields a 
relation for the normalized voltage v(w) in terms of the normalized 


electric field e,(w) and the normalized conduction current i. (w): 


i (w)Le(w)-1] 


v(w) = te ee en)? ; (A.26) 
J(w/w.)D(w,/w.) “[1+(w./wy)] 

Returning to Eq. (A.21) and using 

V(w) = 1, (w)Z(w) ; (A.27) 


where Z(w) is the impedance of the external circuit, we obtain 


Ty (w) = I (w)e(w)-§ (w/w )w Cy(w/wy WE, (w) + 


-J(w/w.)w Ca(w/wy)I,,(w)Z(w) ‘ (A.28) 


Normalizing and introducing the normalized impedance 


Zz(w) = w C4Z(w) (A.29) 
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yields the equation 


: + 7 W a : ie “etd 2 
Veer eco aitaeeny c'ele) iG rE “Wg 


C Le b a 
= 1 (w)e(w)-5(2)0(-8) 2 2e_ (a) 
= Te W My) 7" 7 waa W 
(A.30) 
Elimination of i. (w) from Eqs. (A.20) and (A.30) yields 
| §(0/w,)D(w,/a,)”“L6(w)+(w,/m4) Je, (w) 
Tees O(o)= FF (w/a Jeo lLTF(w, /wy) TT) (A.31) 


Similarly, the elimination of ae from Eqs. (A.20) and (A.30) results 


in an expression for i. (w): 


5 (w/wg)D(0g/og) C145 (w/w,)2 (w) [1+ (Wa/my)] #(u,/ Wg) 3 


i, (v) ¥ 6(w)- (149 (w/w )z(o)T1# (w/w e,(v) 


(A.32) 


Eqs. (A.31) and (A.32) can be formally simplified by introducing 


6(w)+(w./W4) 


v(w) = B(o)-C (aw )z(a)TF(w, JWI (A. 33) 


On substituting Eq. (A.33) into (A.31) and (A.32), we obtain 


ig, (o) = 5(w/u,)(w,/o,) Dv (u)e, (w) (A.34) 


ex a 


and 


i_(w) = 5(w/o,) (w/o) Dl v(w)-Te,(w) (A.35) 
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The ratio of the two currents is 


i, (w)/i (wo) = v(w)/Lv(w)-1] (A.36) 
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APPENDIX-B 
LARGE-SIGNAL DESCRIPTION OF THE READ DIODE 


In a deterministic analysis, the stochastic term in the Read 
diode equation (4.13) is deliberately left out so that this equation 
reduces to Eq. (4.4) which we shall rewrite here, using the symbol 
Tig(t) to denote the deterministic current: 


I 
(2+ (Me) Tr og(t) = (B.1) 


on 
It has been shown in Ref. 20 that, under usual operating conditions, the 
conduction current I g(t) can be expressed by means of the electric field 
component EA in the avalanche region in the following manner: Let fe be 


defined as the ratio of the Bessel functions of order n and zero, respect- 


ively: 
fy = TCD) )/T,(D)) (B.2) 

where 
D, = (m/w kt.) |E.,/E,| = Dle.4| : (B.3) 


Using the expansion 
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cd =] + 5 ) 1 pexP(Jnwt) , (B.4) 
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the normalized conduction current components Cae are given by? 
w 2ee fi 
sa) oa ar! “ su 
het 2(-1) Oe ) ) yen’ pen of fp) Sata (B.5) 
n40 eup=e p 


This relation has been derived for the electric field phase angle chosen 


to be equal to 1/2: 
ei Jjle.y| ; (B.6) 


The normalized electric field ea is the initial value with which 


the large-signal calculations are started. Having obtained the normalized 


conduction current components ae the external current component a at 
the fundamental frequency Wo is obtained from Eq. (A.19): 
Tee aa il) cen ee fo jar le of (B.7) 
ex cl agvaC al cl aeac al : ; 


The normalized diode voltage component Vy° at the fundamental frequency 


We is computed from Eq. (A.26): 


1_,[e(w_)-1] 
gio (0.8) 
J(w,/o.) “DEI tw,/wa] 
In order to sustain oscillations, the external circuit must exhibit the 


normalized impedance (defined by Eq. (A.29)) 
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Z(wo) = we CgZlwg) = welCyVy/Tayy = Cw Cyl o/Va) (Vy/ieyy) = 


= (D/2)(wq/a,) “Viligyy 


If the external circuit is a simple series combination of a 
reactance and resistance, the normalized output voltage component Ve at 


the carrier frequency We is given by 
ye =oVi/VctemR (V1)/V, = R (v1) ; (B.10) 
The total normalized power obtained at frequency We is calculated from 


Pr = PL/I, = 0.5V ia] (B.11) 
To make comparisons with actual diode measurements, the effect of 
parasitic series resistance must be taken into account. The ratio of the 


effective normalized output power and the total normalized power is given 


by 
Poge/Pp = 1 - ro/RQZ(w.) > (B.12) 


where ne is the series resistance normalized with respect to V/w Cy. 
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APPENDIX-C 
DERIVATION OF THE COEFFICIENTS an-2,2 IN THE CONVERSION MATRIX (4.44) 
Using the normalized electric field and current, introduced in 


Appendix-A, Eq. (4.36) may be written in the form 


(Erm strani) (c.1) 


fe) 0. 
= 2 
) Ae, (w nw _) Be ave 


n=-M : 


where the coefficients rN and Bos from Eqs: (4.37) = (4.40), are 


> 
fo) 
1 


= JjweAE, [m(E., = E/E, -dkwt J-m1, (Gae2)) 
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j(mweA/2)E. -(m/2)1 4 


ow 
Oo 
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(nye je | (c.5) 
Cane D ati 
n#0 
Eqs. (C.2) - (C.5) should now be also normalized. First we note 
that 


m(E Eee, alot} : (C.6) 


ao- 


where M is the average (dc) multiplication factor. In analogy with Ref. 8 


we introduce the parameter 
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(6.7) 


e dimension of frequency. Now, using the normalizations and 


parameters introduced in Appendix-A, the coefficients Ne and BP can be 


rewritten as 


o> 
A, = 


= 
oO 
iT] 


n#0 


" 


follows: 
mI { 2 Em . -jw]-1} = 
fe) ml E,kr. 


2 
mI {5 (w/w,) [-(wytjw)]-1} = 
ml {5 (w,/0,)* (w/w) E(uy/o, 45 (w/w )1-1 
mI Ljwk/m + kt wy/m] 


mI Ci(w/w)/D 7 (wy/w)/D] = 


(m1,/D)[(wy/w.) +5 (w/,)] 


(mI,/2)[5(weAE,/ml,)me, -i. J 
(mI /2){§(weAE, kr, /mI,) (wom/u kt) aa Veyn 
(m1, /2)[5(w/w,) (w,/o,) Dea teva] 


= -(mI/2)e,, 


Eq. (A.31), Eq. (C.1) may be written in the form 
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(C.8) 


(C210) 


(C.11) 
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M 


a_(w)e 


a fy (wtnw ) = F(w)/ml : (Gala) 


as 


where a, (w) is determined by means of Eqs. (C.8) - (C.11) and (A.33) - 
(A.36) as follows: 


a,(w) = [AO + BOi(u/w.) (w,/w,) Do(w) (mI) | 


-j(w,/0,) “(w/o,)L(wy/a,)+5 (w/a) 1-1 2 


D”'[(uy/ag)+5(w/ag) 5 (w/og) (ag/ag) Dy(w) = 


j(w/w,)(wq/a,) C(wy/o )+5(w/o,)IEvlw)-1]-1 0. (C13) 


CAP + B° jfornw,)/w 1(w,/u.)“Dv(wtnw.) (ml) = 


a,(w) 
n#0 


(1/2)L5(w/u,)(w4/0,) De. tin(w,/w,) “Dv"(nw Je] + 


Cc 


-(1/2){5L (wtnw )/o T(w,/o,) Dylwtnw ea) = 


wtnw 


: wtnw_)} 
xc "C € eae 
W tw, * 
2 aera saan v(wtnw, )] + nv (nw )-1]} = 
* aly (nw oe 


-(40,/2)14L (wtnw,)/w JE1-v(wtnw,)I/nEy (nw,)-1]) 
(C.14) 


The coefficients an in Eq. (4.44), Chapter IV, are obtained as 


follows: 


a = a, (w=atew) ; | (C.15) 
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with 
p-e =n. (C.16) 


The parameter 


6(w) a (w./W4) 
UREA oy ere sara bo wie oa, )z(w)L1¥ (W/W) (C.17) 


can be somewhat simplified for specific frequencies: 
v(a) = [i(a/w,)z(@)V! (C.18) 


since 6(2) = 1. for 2 << Wes 


= 0 since Z(w>>w)> 
(nw) 


[n|>1 (6.19) 


Finally, from Eq. (A.36) we have 


Wo.) = E-ia/iggl) (C.20) 


The conduction current components iL a and the external current 
joy? needed in the expressions for the coefficients a, (w)» are obtained 


from the large-signal analysis described in Appendix-B. 
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